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Section 1 


SUMMARY 

The primary objective of this multitask program was to continue the development of 
spiral artery cryogenic thermal diode heat pipes. Ethane was the working fluid and stainless 
steel the heat pipe material in all eases. The major tasks included: (1) building a liquid 
blockage (blocking orifice) thermal diode suitable for the HEPP space flight experiment; 

(2) building a liquid trap thermal diode engineering model; (3) retesting the original liquid 
blockage engineering model, and (4) investigating the startup dynamics of artery cryogenic 
thermal diodes. 

The liquid blockage diode for HEPP was U-shaped and had a two diameter stainless 
steel envelope (6.35 mm, evaporator and 9.525 mm, transport/condenser). An electrically 
heated copper block was bolted to the evaporator and a flanged aluminum saddle was brazed 
to the condenser so that a radiator surface could be attached to it. Qualification test re- 
sults were near predictions except for one shortcoming. Start up attempts directly from a 
room temperature cooldown or from an adiabatic condition after a hard thermal dry-out 
were erratic although a few successes were achieved. Successful priming could only be 
effected by first entering a reverse mode (i.e. , flooding the evaporator section), which 
both pressure primed the artery and wetted the evaporator surfaces, followed by a gradual 
condenser cooldown (~ 1° C/Min). A summary of the measured performance for all of the 
heat pipe diodes built under this contract follows. 




Liquid Blockage 

Liquid Trap, 



Eng Model 

HEPP Model 

Eng Model 


(QL) , W-M 

'max 

11.2 

24 

12 

0 

Forward Conductance, W/°C 

4.5 

5 

6 

o 

Reverse Conductance,- W/°C 

0.005 

0.005 

0.037 

e 

Shutdown Energy, W-hr 

0.02 to 
0,10 

0.10 

0.32 


Since the startup difficulty could only be overcome by undergoing a preconditioning 
sequence, this might limit the possible application of artery thermal diodes. An investiga- 
tion was therefore undertaken to study and possibly resolve the problem. Two reasons 
were suspected: artery priming difficulties and/or surface wetting problems. An analytical 
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investigation of artery priming determined that artery permeability is the important 

criterion for how easily an artery primes. For a given capillary pore size, the higher the 
wick permeability the larger the evaporator heat load at which hill artery priming will occur. 

An experimental investigation was also conducted into the wetting characteristics of 
ethane/stainless steel systems using a specially constructed. chamber that permitted insitu 
‘observations. Results were qualitative but revealed that there are no wetting problems 
with the ethane/stainless steel system. Ethane exhibited the same predictable wetting be- 
havior as Freon 113 and Acetone which were used as comparison standards. 
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Section. 2 


INTRODUCTION 

The application of heat pipe thermal diodes to the cryogenic temperature range would 
greatly expand the mission capabilities of the ever-increasing number of spacecraft that 
require low temperature cooling for their payloads. Instead of relying upon solid cryogens 
or refrigerated helium loops, which are either heavy or unreliable, the thermal diode would 
permit the use of a passive space radiator system for almost any vehicle orientation. -When 
the radiator surface becomes cold enough to serve as an acceptable heat sink (e.g. , views 
deep space), the diode heat pipe would operate as a high thermal conductance link between 
the radiator and the heat .source. This allows the payload to be cooled by radiation. If the 
radiator surface becomes, warmer than the heat source, as when illuminated by sunlight, 
the diode action would break the thermal connection thereby isolating the -payload. During 
this time, the payload would either remain cold due to its own thermal inertia or it could be 
coupled to an alternate heat sink such as a solid cryogen. In any event, the benefits of a 
viable heat pipe diode thermal control system are obvious - less overall system weight 
and/or longer mission life. 

The thermal diode concepts that are considered the most practical for hardware 
implementation because they result in the fastest, most positive shutoff, are the liquid trap 
and liquid blockage configurations. 

In the liquid trap design, a non-communicating reservoir is located at the evaporator 
end of the heat pipe and also is in thermal contact with the heat source. During normal 
mode operation, when the evaporator and reservoir are warmer than the condenser section, 
the reservoir is empty because any fluid in it has vaporized. The fluid cannot return since 
there is no wick communication with the evaporator. When the thermal gradient reverses 
(condenser warmer than evaporator) heat pipe operation is terminated because all of the 
working fluid eventually condenses and stays trapped in the reservoir. Thus, the heat pipe 
wickbecomes starved of fluid, and cannot operate. Because the liquid trap diode mechanism 
relies on complete removal of all working fluid, the shutoff time is directly related to fluid 
inventory. It is usually better suited to those heat pipe wicking systems that require less 
fluid (i.e., axially grooved heat pipes as opposed to arteries), since the shutoff energy 
requirements are lower and shutoff time quicker. 
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The liquid blockage concept also uses a non-communicating fluid reservoir but with 
two important differences. The’ reservoir is located at the condenser end of the heat pipe 
and it is completely filled with excess fluid during the normal operating mode. Shutoff 
is effected when, due to the reversed temperature gradient, the excess fluid evaporates 
from the reservoir* and condenses in the vapor space of the evaporator section. The 
. reservoir is sized so that it holds enough, fluid to completely choke the evaporator vapor 
space, which interrupts the heat pipe action. 

The conventional liquid blockage design needs a vapor space small enough to self- 
■ prime during 1-g testing. But this restriction also results in decreased transport capacity 
during normal operation due to excessive vapor pressure losses. A unique concept which 
permits blockage without requiring a restricted vapor space is the blocking orifice design 
(References 1 and 2). In this concept, a solid plate with one small opening (orifice) adjacent 
to the wall, is located inside the heat pipe at or before the entrance to the evaporator sec- 
tion. During normal operation, a small vapor pressure drop is experienced at the orifice. 
During the reverse mode blockage is established by die creation of a liquid meniscus across 
the orifice, which retains (without draining) all of the excess liquid which has condensed in 
the vapor space of the evaporator. 

As previously mentioned, a very common aerospace application of thermal diodes in- 
volves heat rejection from a source package to a radiator during normal operation, with 
shutoff when the radiator is exposed to a hot environment. The conventional liquid blockage 
’ diode concept or liquid blockage using an internal orifice are attractive shutoff techniques * 
in this case for the following reasons: 

■ ® The liquid reservor can he easily heated by the hot environment to expell liquid 
in the shutoff mode and it can he cooled by radiation to a cold environment to 
promote liquid retention in the normal mode 

9 Generally, the evaporator is short compared with the condenser and total pipe 
length, minimizing blocking fluid requirements 

© The liquid reservoir volume may be substantially reduced or eliminated by the 
extent to which the blocked length liquid requirement is satisfied by expansion of 
the normal fluid inventory during the hotter reverse mode condition. 
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The liquid trap concept for blockage in this case is less attractive because of 
packaging considerations. For reverse-mode operation, the reservoir to collect the normal 
inventory of working fluid must become an integral part of the heat source package. Thus, 
not only must there be a thermal interface between the heat pipe evaporator and source but 
also between the source and the heat pipe reservoir, which complicates the source package 
assembly. 

This report describes four major tasks which were aimed at advancing the development 
of spiral-artery heat pipe thermal diodes. 

1. Build and test a liquid blockage (blocking orifice) cryogenic thermal diode for the . 
HEPP experiment 

2. Fabricate a liquid trap cryogenic thermal diode engineering model for testing at 
NASA- ARC 

3. Thermal vacuum test the original liquid blockage (orifice) cryogenic thermal diode, 
engineering model 

4. Investigate (by analysis and experiment) the start-up dynamics of spiral artery 
cryogenic thermal diode heat pipes. 
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Section 3 


FABRICATION AND TESTING OF A LIQUID BLOCKAGE (BLOCKING ORIFICE) 

CRYOGENIC THERMAL DIODE FOR THE HEPP SPACEFLIGHT EXPERIMENT 

(MODS 8 AND 10) 

3.1 THE EXPERIMENT 

The main obj ective of the HEPP spaceflight experiment is to measure long term 
zero-g operation of selective cryogenic heat pipes. The experiment incorporates two types 
of heat pipes : a conventional axially grooved (rectangular groove profile) aluminum iso- 
thermalizer, and a stainless steel thermal diode. Both heat pipes use ethane as the cryo- 
genic working fluid and have a nominal operating temperature around 200 K (-100° F). 

The basic experiment shown in Figure 3-1, is mounted in a tray and consists of a 
shielded space radiator, the two heat pipes, and a phase change material (PCM) canister 
which contains heptane (melt point = 182 K). For packaging purposes the heat pipes are 
bent in a U-shape, the condenser sections are mounted to the radiator and the evaporators 
are positioned inboard and attached to heater blocks. The PCM canister, which has a 
nominal 30' watt-hour capacity, is also mounted to the radiator and serves as the heat sink 
for the high power tests, since the dump capability of the radiator is small («s 3 watts) and 
suited only for the long duration, low power tests. The radiator also serves to resolidify 
the heptane so that it can be used again. 

Two kinds of thermal diode concepts were considered: a liquid trap design (made with 
axially grooved forged stainless steel tubing), and a liquid blockage design which used a 
spiral artery wick and the blocking orifice technique. This section describes the develop- 
ment of the liquid blockage thermal diode concept. 

3.2 DESIGN AND FABRICATION 

Three candidate envelope configurations were considered for the cryodiode heat pipe: 
(1) the 1/4 inch (6. 35 mm) OD tube used for the engineering model (see Section 5), (2) 
a larger 3/8 inch (0,53 mm) OD tube, and (3) a hybrid two-diameter design that has a 1/4 
inch (6.35 mm) OD evaporator /blocked-transport section and a 3/8 (9.53 mm) inch OD con- 
denser/uhblocked-transport section. The baselined diode concept for each one was the 
blocking orifice-liquid blockage technique using ethane as the working fluid, which was 
successfully demonstrated in the engineering model. 
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Fig. 3-1 HEPP 
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After running various parametric analyses which altered the artery diameter, it was 
decided that the most predictable design would be one that used the same artery developed 
for the engineering model. In this way existing tooling and proven fabrication techniques 
could be directly utilized. Performance predictions were made for each candidate configura- 
tion using the same 0.117 inch (2.97 mm) OD artery. Calculated throughput for the 3/8 
(9. 53 mm ) OD tube and the hybrid 2-diameter designs were nearly equal (9. 5 watts @ 190 K, 
o-g) and about 10% higher thanthe 1/4 inch (6.35 mm) OD configuration. 

However, besides throughput other things must be considered. A prime disadvantage' 
of the uniform l/4 inch (6.35 mm) OD design is the relatively narrow vapor annulus (0. 0385 
inch (.978 mm)) that is kept throughout the unblocked transport section and the condenser, " 
when it is only needed in the evaporator section. The narrow vapor space in the evaporator 
is required since it permits testing in one-g and also minimizes the excess liquid and 
reservoir volume requirements. In addition to creating assembly problems, maintaining 
a narrow vapor space in the condenser increases the likelihood that the vapor annulus will 
self-prime and interfere with normal forward mode performance. 

Using the larger 3/8 inch (9.53 mm) OD tube eliminates the self-priming problem but 
introduces another disadvantage. It requires a larger fluid reservoir (3.9 times larger) to 
hold the excess fluid needed to block the evaporator vapor space during the reverse mode. 

The hybrid two-diameter configuration offers a reasonable compromise. It has a 
slightly higher throughput over the constant 1/4 inch (6.35 mm ) OD envelope, while main- 
taining the relatively small reservoir size and also providing the advantage of a larger vapor 
space through the condenser. A comparison of the three concepts with respect to reverse 
mode considerations (reservoir size, shut-down losses) is presented below for a 190 K 
evaporator and a 210 K condenser/reservoir. 

Reverse Mode Comparison (Ethane Fluid) 


Design Concept 

HP Charge 
(grams) 

Charge 

(grams) 

Reservoir 

Volume 

(in 3 ) 

(1) 

Length 

(in.) 

Shut down Energy (2) 
during 1 hour 
(watt - hr) 

1) 

1/4 OD 

3,78 

1.13 

0.127 

2.0 

0.19 


(6.35 mm) 



(2.08 ec) 

(50.8 mm) 

2) 

3/8 OD 

5.3 

4.42 

0.496 

7.6 

0.65 


(9. 53 mm) 


w 

(8.13 cc) 

(193.0 mm) 

3) 

1/4-3/8 OD 

5.0 

1.13 

0. 127 

2.0 

0.19 


(6. 35-9. 53 mm) 


(2.08 cc) 

(50. 8 mm) 
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(1) Reflects engineering model-type design (6.25 in. OD) (158.75 mm) 

(2) Includes conduction and transition energy losses. 

In addition.to the above, the pipe pressure at storage conditions (120° F) (322 K) is 
much lower for the hybrid design than the others; 886 psia vs. 1449 psia (6108 KPa vs 
9989 KPa) for .1/4 OD and 1036 psia (7142 KPa) for 3/8 OD. 

Since the hybrid 1/4-3/8 OD design offers the best combination of forward mode and 
reverse mode performance, it was selected as the baseline configuration for the liquid 
blockage diode heat pipe. An overall layout of the basic cryodiode configuration is shown 
in Figure 3-2. Performance predictions as a function of operating temperature are given 
in Figure 3-3. 

Early in the program several hardware samples were fabricated to verify proposed 
assembly techniques. 

o A sample artery was made and tested by performing lift tests before and after -* 
several bending cycles. The lift tests verified no degradation in pore size due 
- to bending 

9 3/8 to 3/8 (9. 53 mm to 9.53 mm ) OD butt weld samples made with an orbital 
welder showed good weld penetration and negligible drop-through 

e Swaging and spinning operations for the transition sections were .also successfully 
demonstrated. Samples were made which join l/4 inch (6.35 mm) OD (evaporator) 
to 3/8 (9. 53 mm ) OD (condenser) tubes and 5/8 (15. 88 mm) OD (reservoir) to 
3/8 (9. 53 mm) OD (condenser) tubes. Both samples were butt-welded (orbital) 
together with good results 

3.3 RESERVOIR DESIGN CONSIDERATIONS 

A non- communicating reservoir is usually preferable to one that can communicate 
with the main artery. Although the latter has the advantage of being able to replenish the 
main artery in the event of a fluid deficiency, it also has other counterbalancing 
disadvantages. 

• If the artery wick is allowed to extend beyond the end of the condenser, it causes 
a backpull on the condenser meniscus which decreases the net available capillary 
head and results in a lower transport capacity. Since cryogenic heat pipes are 
typically low capacity pipes to begin with, this further decrease in capability is 
undesirable. 
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Fig. 3-3 Cryodiode Forward Mode Predictions 
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• A communicating reservoir requires extra wicking which in turn occupies volume 
that must be compensated for by increasing the reservoir size significantly. This- 
is especially time with our reservoir design that has individually separated 
capillary channels and a highly conductive aluminum- matrix to accelerate evapora- 
tion in the reverse mode. 

® The extra wicking in a communicating reservoir increases the required fluid 
charge, which in turn increases the storage pressure.. 

Temperature control requirements for the reservoir of a liquid blockage diode are 
relatively simple: the reservoir should rim slightly' (a fraction of a degree K) below vapor 
temperature in the normal mode of operation, and slightly above vapor temperature in the 
reverse mode of operation. Both requirements can -be met by thermally coupling the reservoir 
to the condenser boundary conditions. 

As an example for illustrative purposes, consider a spacecraft application in which a- 
diode couples a detector to a radiator which is subj ect to a varying space environment. The 
diode condenser is attached to the radiator. The excess liquid reservoir is at the condenser 
end of the diode, and has a small portion of the radiator surface dedicated for reservoir 

t ' ’ * 

thermal control. During normal mode operation, the reservoir will tend to reject heat 
from the dedicated radiator surface. The axial solid conductance between the condenser 
and reservoir is relatively small. The reservoir can, however, obtain heat by condensation 
if only partially filled with liquid. When the reservoir internal passages are filled with 
liquid, and a thin liquid layer, covers the exposed face of the reservoir internal aluminum 
block, the thermal resistance between the vapor and the reservoir will receive only a 
small amount of heat by condensation, and the reservoir temperature will fall below the 
main radiator temperature, to a value approaching the adiabatic surface temperature. 

The preceding discussion involved three assemptions: 

1) negligible heat leakage from the spacecraft to the- reservoir 

2) equal absorbed heat flux on the reservoir radiator and main radiator 

3) negligible variation in diode boundary conditions during normal mode operation. 

Heat leakage must be minimized.between spacecraft and diode and spacecraft and radiator 
for any reasonable cryogenic temperature. Insulation requirements between reservoir and 
spacecraft will be similar to those for other parts of the diode and radiator, and well 
within the capabilities of multilayer insulation blankets. Equality of absorbed, heat flux 
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should occur naturally, since both radiators will usually he sharing the same surface. The 
rate of variation in diode boundary conditions is a function of orbit and surface orientation. 

Considering the same example under reverse mode operation, when the absorbed 
flux drives both radiator surfaces above diode vapor temperature, evaporation occurs in 
both the reservoir and condenser regions and condensation occurs at the detector end of 
the pipe (evaporator). Excess liquid obtained from the reservoir will accumulate at the 
detector end of the diode, effectively shutting off the diode. As the diode shuts off, the 
temperature of the vapor, condenser, and reservoir all rise. To keep the reservoir above 
vapor temperature, the inequality: 

/ Q abs \ 

\ MC p ) res 

must be satisfied. Reservoir temperature will then rise faster than condenser temperature, 
but will stay close to vapor temperature so long as liquid is available in the reservoir for 
evaporation. 

For applications in which (he radiator environmental flux varies rapidly, the transient 
cool down of the rdservoir following a shutdown must be at least as fast as the condenser to 
permit the reservoir to fill with liquid when the diode is ready to turn on. This implies 

f— ) 

\ MCp J res 

Where A R for the reservoir is that portion of the radiator which is segmented and 
dedicated to reservoir cooling, and A^^ for the condenser is the main radiator area. 

The MC product includes the entire mass which must be cooled by heat rejection to space 
(i. e. for the reservoir - the reservoir case and internal metal, attachment saddles, and 
reservoir radiator; for the condenser - the condenser including wick and fluid, the attach- 
ment saddles and the main radiator). 

It may be noted that the temperature eontrol requirements for a liquid blockage diode 
are relatively easily achieved compared with those for the gas reservoir variable conductance 
heat pipe (VCHP). The gas reservoir (or gas reservoir cold trap) must in general be main- 
tained at a temperature significantly below condenser temperature to permit gas storage at 
a reasonable gas partial pressure in the reservoir. Hence VCHP's are subject to changes 
in set point as a result of coating degradation or anticipated thermal leakage. The liquid 
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reservoir, by contrast, runs close to condenser temperature. Thermal isolation of the 
reservoir from the spacecraft can be handled together with the rest of the diode and radiator, 
since requirements are similar, and the coating degradation constraint is no more severe 
for the liquid reservoir radiator then it is for the main radiator. 

The ability of the reservoir to function properly with small temperature differences 
has been confirmed experimentally. A l/4 in. (6.35 mm) liquid blockage diode' with ethane 
working fluid was operated continuously for 5 hours at 10 watts evaporator power -at 1 cm - 
adverse tilt and a temperature of -100° F (200 K) with the reservoir at a temperature only 
0.02 to 0.04° F below vapor temperature. 

‘3.4 ORIFICE ORIENTATION 

The orientation of the orifice has an effect on the blocking capability of a given size 
opening under 1-g test conditions. Two orifice orientations are of interest: 12 o'clock and 
6 o’clock. The maximum tilt position at which complete blockage can still be realized is 
obtained by equating the orifice capillary pressure to the hydrostatic head that it must 
support. 

3.4.1 12 O' Clock Orifice 



P^ = Pressure at lowest point in condenser liquid (maximum pressure point, 
° £Lat menescus) 

P v = Vapor pressure 



Orifice Capillary Pressure 

a 20r 

AP = = = (P -P ) 

_ho ho v i 12 

2 
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These results have been plotted in Figure 3-4 as a function of both heat pipe temperature 
and adverse tilt. As seen, the analysis predicts a more tilt sensitive orifice at the 12 o’clock 
position, which is not intuitively obvious since it is at the top of the pipe. However, this 
paradox can be explained by studying the hydrostatic balances developed and realizing that 
the blocking meniscus must sustain a larger vapor-liquid AP at the 12 o’clock position 
since the liquid pressure is lower than the 6 o’clock point, while the vapor pressure re- 
mains the same. For adverse tilts, complete blocking (up to the orifice plate) at the 12 
o’clock position can only occur at the colder temperatures ( ~ 150 K). However, blocking 
throughout the 4-inch (101.6 mm) evaporator section only, will take place over the entire 
temperature range from 150 K to 210 K. The 6 o’clock position can sustain complete 
blocking up to the orifice over the 'entire temperature range for adverse tilts below about . 10 
inches (2.54 mm). 

Although the preceding discussion was based on the analytical predictions, test data 
have been observed that support the results. An evaluation of the test data which was ob- 
served during the checkout runs follows, 

e At the 12 o’clock orifice position and .10 inch (2.54 mm) adverse tilt only the 4- 
ineh (101. 6 mm) evaporator section was successfully blocked at an evaporator 
temperature of 200 K 

• At the 6 o’clock position, complete blocking (up to the orifice plate) was observed 
for level, .05 (1.27 mm) and .10 inch (2.54 mm) adverse tilts at a 200 K evaporator 

• With the orifice still at 6 o’clock, when the tilt was further increased to .15 inches, 

, (3. 81 mm) blocking at the orifice could not be achieved. Evaporator temperature 

was also 200 K. .... ’ 

3.5 THERMAL MODEL 

The simplified thermal nodal model of the liquid blockage diode heat pipe experiment 
which was developed is presented in Appendix A, including the AGTAP subroutine code for 
the diode mechanism. In addition to altering forward and reverse mode parameters in 
accordance with the calculated temperature gradients, the program logic includes a PCM 
routine which models the phase change phenomenon. It also has generalized inputs to permit 
rapid modification and evaluation of radiator dimensions and location relative to diode attach- 
ment points. Two options exist for .controlling the inputs to the PCM and reservoir heaters. 
One uses an on/off duty cycle routine which permits separate changes in both heater wattage 
and duty cycle. The other heater control option uses a fixed temperature set-point and 
tolerance band. When the control temperature is above the set point, the heater is off; when 
it is below it,, the heater is turned on. 


3-11 



3-12 




ADVERSE 


FAVORABLE 


{alp) ACTUAL 


ASSUMES BLOCKED EVAPORATOR TEMP 36°F (20 K) 
COLDER THAN VAPOR TEMP 


/5@T e ,ct@T e + 36 0 F 



-190 -180 -170 -160 -150 -140 -130 -120 -110 -100 -90 -80 -70 

•BLOCKED EVAPORATOR TEMP, °F 


0787-004 


TILT, IN 


210 K 





Computer runs were made which verified that during the forward mode the reservoir 
temperature will always be lower than the condenser (and vapor) temperatures by about 1 K, 
thereby ensuring a filled reservoir and clear' condenser vapor space. The predicted thermal 
profile after 4 hours with an evaporator input of 10 watts is given in Table 3-1, Runs were 
also made to demonstrate successful reverse mode transition and temperature separation 
between condenser and reservoir during the following cooldown. Results are given in 
Figure 3-5 for the case of a 40 watt PCM input (40% duty cycle) and a steady 2 watt reser- 
voir input. 

The effect of reservoir heat input on blocking response was also investigated at two 
extreme PCM heat inputs (5,6 and 40 watts). The results are given in Table 3-2 and indicate 
a minimum reservoir heater requirement of at least 2 watts to achieve complete blockage 
in a reasonable time (~10 minutes). . 

3.6 HARDWARE DESCRIPTION 

The HEPP diode heat pipe utilizes the blocking orifice technique to effect liquid 
blockage of the evaporator vapor space during reverse mode operation. It is made- from 
‘ stainless steel tubing (28 mil wall) and has a l/4 inch (6.35 mm) OD evaporator /blocked 
transport section, a 3/8 inch (9.53-mm) OD condenser/transport section, and a 5/8 inch 
(15. 88 mm) OD excess fluid reservoir. The different sections are joined together by 
: swaging the appropriate tubes to achieve a uniform interface OD and then butt welding. 

The pipe contains a spiral artery wicking system similar to that.used in the engineering-model 

t 

(see Section 5). and uses purified ethane as its working fluid. 

As shown in Figure .3-2, the diode is bent into a U-shape, One leg of the U contains - 
the evaporator, which is attached to a cylindrical copper heater block; the opposite leg con- 
tains a brazed aluminum saddle which is attached to a radiator surface. The excess fluid 

' reservoir is located at the end of the condenser and also contains a brazed aluminum saddle 
2 2 - 

to which a 9 in (58 cm ) radiator surface is attached. The overall dimensions of the U- 
shape are approximately 14 inches (0.36 m) across the base and 22 inches (0.56) along each 
leg. A detailed weight breakdown is contained in Table 3,3. 

The finished diode heat pipe assembly is shown in Figure 3-6 after the saddles were 
.-brazed to the condenser and reservoir, and the mounting tabs welded onto the transport 
section.. A closeup view. of the. reservoir assembly which shows the clips for attaching the 
PRT sensors is shown in Figure 3-7. Closeups of the brazed condenser assembly and the 
transport-to-evaporator transition section are shown in Figures 3-8 and 3-9. Note the 


3-13 



TEMP, K 


Table 3-1 Normal Mode Experiment Temperature at 4 hours,- = 10-Watts 


NODE 


TEMP (K) 

1 

. EVAPORATOR BLOCK 

193.6 

2 

EVAPORATOR WALL 

192.7 

3 

HP VAPOR 

192.0 

4 

CONDENSER WALL 

191.9 

5 

RESERVOIR 

191.1 

6 

RESERVOIR RADIATOR 

1909 

7 

MAIN RADIATOR (W/PCM) 

183.3 

8 

MAIN RADIATOR 

191.6 

g 

PCM CANISTER 

182.6 

0787-035 


* 


. -3W0N RES 


■CONDENSER 


' CM- .5 HR 40W 

> .5 HR 0 - 

DUTY CYCLE : PCM = 40% 
RES = 100% 


.RESERVOIR AREA - 9 IN. 2 (58 cm 2 ) 


RESERVOIR' 


0787-005 


- TIME, HOURS 


Fig. 3-5 Reverse Mode and Cool Down Predictions for HEPP (Cryodiode) 
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Table 3-2 Reverse Mode Performance Predictions 


INITIAL PCM 
TEMPERATURE 


HEAT INPUT TO 
RESERVOIR (WATTS) 


HEAT INPUT TO 
PCM (WATTS) 


TIME FOR BLOCKAGE 
TO OCCUR (MIN) 


BLOCKED 
LENGTH (INCHES) 


-100 

1.0 

5.6 

15.6 

-100 

2.0 

5.6 

62 

-100 

5.0 

5.6 

2.4 

-100 

1.0 

40.0 

4.6 

-100 

2.0 

40.0 

9.8 

-100 

5.0 

40.0 

2.6 

-162 

1.0 

40.0 

31.7 

-162 

2.0 

40.0 

8.9 

-162 

5.0 

40.0 

2.6 


PARTIAL 

BLOCKAGE 


PARTIAL 

BLOCKAGE 


8.77 

8.80 

8.85 

-1.04 

8.88 

8.86 

-4.00 

8.18 

8.16 


0787-036 


Table 3-3 Detailed Weight Breakdown, HEPP Liquid Blockage Cryodiode 


COMPONENT 

MATERIAL 

WEIGHT (GRAMS) 

1) HP TUBE ENVELOPE 

304 STAIN LESS STEEL 

36.7 

2) SPIRAL ARTERY 

• 304 STAINLESS STEEL 

32. 

3) RESERVOIR INSERT 

ALUMINUM-6061 

25.4 

4) WORKING FLUID 

ETHANE 

10.2 

5) EVAPORATOR BLOCK 

COPPER 102 

998.8 

6) SADDLE ASSY 

ALUMINUM-6061 

272.4 

7) RESERVOIR RADIATOR ' 

ALUMINUM-6061 

45.4 

8) INSULATOR SHIM 

KAPTON 

5. 

0787-037 


« 



Fig. 3-6 Spiral Artery Liquid Blockage Cryogenic Thermal Diode for HEPP 


3-15 













3-16 






0787-009 


Fig. 3-8 Brazed Condenser Assembly 
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Fig. 3-9 Two Diameter Transition Section 


3-17 














3-18 










PRT clip and mounting pads. Figure 3-10 shows the final diode heat pipe assembly mounted 
to its vibration test platform. At this point, it had been charged with ethane, welded closed 
and thermal vacuum tested. 

3.7 QUALIFICATION TESTING 

Objective. The main obj eetive of the qualification test was to demonstrate a diode 
heat pipe that could be effectively used as part of the HEPP experiment. Consistent with 
this objective, the pipe requirements were: ,a minimum forward mode throughput of 750 

i , 

watt-cm (0— g equivalent) and good thermal conductance, at least 5 watts/ 0 C. It should also 
show a rapid transition to the shut-off operating mode in response to heat inputs at the con- 
denser and reservoir and quick recovery to forward mode operation. . 

As a secondary objective the effect of orifice orientation on 1-g performance in both 
the forward and reverse operating modes was also evaluated. Two orifice orientations 
were considered, (a) orifice opening at the top of the pipe (12 o'clock-position) and (b) 
orifice opening at the bottom (6 o'clock position). As shown in Section 3.4, the 12 o’clock 
orientation is the more conservative one as far as maintaining a blocked evaporator, since 
it must sustain a larger delta-P, - But even the 6 o’clock orientation imposes a greater demand 
than the actual zero-g flight requirement. As far as influencing the reverse-to -forward mode 
transition time, the 6 o’clock position.would be most like the zero-g situation since the favor- 
able draining effect would be less than that experience with a 12 o'clock orifice opening. 

Test Set-up 

The cryodiode heat. pipe qualification unit was tested in a 2 ft (0.61 m) diameter x 4 ft 
(1.22 m) long thermal vacuum chamber located in Grumman's Thermal Laboratory. A 
cross-section of the installation is shown in Figure 3-11. The U-shaped heat pipe was 
attached to a rigid mounting platform through low-k thermal isolators, and the entire set-up 
was isolated from the room ambient by an LN^ cold wall and a 10 torr vacuum. 

Heat was supplied to the pipe through' a 1.5 inch (38.1m) diameter copper heater 
block attached to the evaporator section. At the condenser, heat was removed through an 
aluminum radiator fin (2.6 ft ) (0,24 m ) which was bolted to the condenser saddle. The 
radiator fin contained a tube through which LN 0 could flow to accelerate the system cooldown 
from room temperature. Once near the desired operating temperature, the LN„ flow was 
stopped and the pure radiation coupling to the cold wall used to reject the heat. Strip heaters, 
attached to the back of the radiator fin and insulated with MLI, were used to adjust the heat 
pipe vapor temperature as required. 
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Fig. 3*1 1 Cryodiode T/V. Test Installation (Qualification Unit) 
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The entire mounting platform was balanced between a center support, which ran the 
length of the platform, and an eccentric cam which was positioned .off-center near the 
chamber and door. The cam was attached to a rotary feedthru so that it could be adjusted 
under test. The amount of angular rotation of the plate was measured by a calibrated 
scale which was hard mounted to the plate and which could be sited through the chamber 
door viewing port by using a surveyor’s transit. 

Since the heat pipe was tested with two different orifice orientations (12 o'clock and 
6 o’clock) it was necessary to flipthe heat pipe 180 degrees and remount it to the base 
plate. When this was done, the radiator fin was also reinstalled on the condenser flange 
to' conform to the cold wall shape. 

, Instrumentation ■ 

The qualification unit was instrumented with the 11 platinum resistence thermometers 
baselined for the HEPP experiment and 18 additional copper-eonstantan thermocouples. 

The thermocouples were -used - to check the calibration of the PRTs and provide more de- 
tailed thermal profiles. The PRTs were mechanically fastened to the pipe by suitable 
designed clamps whereas the thermocouples were spot welded to the steel .tubing whenever 
possible. Detailed instrumentation information is given in Figure 3-12.' After instrumenta- 
tion, the entire heat pipe except for the 'condenser and reservoir radiator interfaces, was 
insulated with 25 layers of aluminized mylar (MLI). 

Before starting the formal test procedure, some preliminary test points were run to 
determine if the heat pipe contained the proper charge of working fluid. For these prelimin- 

■ V ' 

ary tests, the charge valve was removed from the pipe by pinching off and welding the 
secondary charge tube, which separated the primary charge tube and the charge valve. After 
the preliminary test data had verified a properly filled'pipe, the primary charge tube- was 
pinched off and welded closed. 

Test Results 

The test program consisted of a detailed performance evaluation before vibration 
testing, followed by selective forward 'and reverse mode spot checks afterward. The 
majority of the tests were run with the orifice at the 12 o'clock position; only during the 
pre-vibration sequence were a few tests run with the orifice at 6 o’clock. 

The forward mode, steady state tests determined the maximum heat pipe throughput 
for a combination of four nominal pipe temperatures (150, 170, 190, and 210 K) and four 
adverse tilt settings (0. 25, 1. 0, 2. 0 and 4. 0 cm). The transient tests, which were run at 



s w 
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Cryodiode Heat Pipe 






the lowest (0. 25 cm) tilt, monitored the reverse mode response for selective reservoir heat 
inputs (1, 2 and 3 watts) at the two extreme temperature conditions (150 and 210 K). 

Forward Mode 

During the forward mode (steady-state) tests the maximum heat load was limited by 
the dump capability of the condenser radiator. The radiator capacity calculated for the 
various test temperatures is given below: 

Test Temp (K) Radiator Capacity (Watts) 

150 4.9 

170 8. 5 

190 12.4 

■ 210 18.5 

At all of the test temperatures, the evaporator heat load could be increased to the capacity 
limit of the radiator, beyond that the heat pipe temperature would increase to accommodate 
the load. Even at adverse tilts up to 4 cm (the limit of the test setup), the dryout limits 
could not be reached. A summary of the steady state performance, in the form of evaporator 
• heat load versus adverse tilt, is presented in Figure 3-13. At the maximum tilt (4 cm) the 
measured heat loads were: 4 watts at 150 K, 7.2 watts at 170 K, 10.5 watts at 190 K and 
17. 3 watts at 210 K. Using a static wicking height of 12 cm (corresponding to lift test data) 
a maximum zero tilt heat load of 25 watts was estimated at 210 K. Temperature profiles 
for representative test points are shown in Figure 3-14. Data are presented for both 
orifice positions, for pre-and post-vibration test points, and for adverse tilt extremes. 

As seen by comparing TP154 and TP104, the orifice position has no effect on the forward 
mode performance. The post-vibration thermal test data (TP101V, TP104V) also show 
little change. The average value. of forward mode thermal conductance is 7.4 watts/ 0 C 
at the 200 K level, and 2.4 watts/ 0 C at the 150 K level. 

Reverse Mode 

Temperature profiles at the beginning and end of -reverse mode tests for the 6 o’clock 
and 12 o’clock orifice positions are shown in Figure 3-15. These fully established profiles 
are shown after 15 and 13 minutes respectively; however, the actual time at which blockage 
was first established was about half that shown. This data is significant because it verifies 
the different temperature profiles predicted for the two orifice locations. The 12 o’clock 
orifice was previously shown to be more sensitive for sustaining a fully blocked vapor space 
up to the orifice location. This is seen in the TP205 data since the blocked section extends 
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, Fig. 3-15 HERP Liquid Block Cryodiode, Reverse Mode Test Results 
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only through the evaporator section (about 10 cm). However, at the 6 o'clock position 
(TP252) the blocked section extends up to the orifice (about 20 cm) and the evaporator is 
more effectively insulated from the unblocked vapor space, which is filled with warmer 
vapor. 


Shutdown energy was estimated by using the method of Reference 4. The shutdown 
time is taken as the interval from the moment the reservoir begins to rise -in temperature 
above the transport section to the time the blocked transport section rate of temperature 
increase drops sharply. The energy absorbed (shutdown energy) during transition is de- 
fined as 

. Q - (mc p ) (T, - T.) - Q 0 T 

where: mc^ = thermal mass capacitance of the diode evaporator and detector block 
Tj; = temperature of the detector block at complete shutdown 
=. temperature of the detector block at the start of reversal 
* Q = detector block heater power - - 

G 

T = shutdown time - • 


This relationship is a conservative estimate since it' includes the energy absorbed by 
the evaporator block due to the surrounding ambient. ^ = . 


’ Being conservative and using the data shown (Figure 3-15) for TP252, ,the shutdown 
energy is calculated to be 0. 098 W-hr at the maximum. ’ , 

The comparison between the platinum resistance thermometers (PRTs) and the 
corresponding copper-constantan thermocouples is shown in Figure 3-16. The PRT calibra- 
tion curve is given in Figure 3-17. As seen, the agreement- between the measurements is 
poor ranging from 0.1 to 5° C difference depending upomthe sensors*. The thermocouple 
measurements are considered the more accurate, since they were spot welded to the heat 
pipe tube. The PRT measurements oh the other hand are influenced by the uncertainty 
associated with the interface thermal resistence between. the sensor and the tube. They use 
mec h anical clips and an interface grease. 

System Testing 


Thermal vacuum testing'of the- integrated HEPP 1 cryogenic heat pipe- experiment was 
conducted at the Johns Hopkins Applied Physics Laboratory (APL) under the direct super- 
vision of ITE, Inc. , the system integrator. Two thermal tests were run; one before and 
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Fig. 3-16 Temperature Comparison PRT vs Copper Constantan Thermocouples 
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one after vibration testing the entire assembly. During each sequence, the diode was tested 
to determine the following: 

• Startup from ambient 

• Long term, low load forward mode performance 

• Maximum heat transport capacity 

• Reverse (shut-off) mode behavior. 

Start Up From Ambient 

Initial start up from ambient was successfully demonstrated during both pre- and 
post-vibration thermal tests by the convergence of evaporator and condenser temperatures 
during the prolonged cooldown. A low evaporator heat load was applied at all times; 0. 5 
watt during the pre-vibration test and 1. 0 watt during the post-vibration test. Figure 3-18 
shows the startup’ data for the post-vibration test sequence. Initially, a large temperature 
difference developed between the evaporator and condenser, but the heat pipe eventually 
recovered. This is indicated by the convergence of the evaporator and condenser tempera- 
tures at -40° C (233 K). ' 

Long-Term, Low Load Performance 

The cryodiode heat pipe was successfully operated at a 2 watt load for several hours. 
Figure 3-19 shows a typical temperature profile taken during pre-vibration tests. 

Maximum Heat Transport Capacity 

The maximum heat transport capacity tests were run at the PCM melting point (182 K) 
by incrementing the evaporator heat load until a dryout was recorded, as evidenced by a 
runaway evaporator temperature. The maximum heat transport capacities recorded during 
the tests are summarized below and reflect operating temperatures from 190 to 200 K. 

Summary of Maximum Measured Transport Capacity 


Test 

Adverse Tilt 
cm (in.) 

Operational 

Q , Watts 

max 

Pre-Vibration 

.19 (.075) 

20 

Post- Vibration 

.19 (0.75) 

24 


.51 (.20) 

22 
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Fig. 3-18 Diode Start-up Behavior During Post-Vibration Testing (1 W on Evaporator) 
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Fig, 3-19 Long Term, Low Load Test Data 
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These data indicate that there was no performance degradation due to vibration and 
that the artery diode heat pipe is very insensitive to adverse tilt. A typical temperature 
profile is shown in Figure 3-20. By comparison, it is of interest to note that the forward 
mode throughput measured for the axially grooved transport heat pipe degraded from 22 
watts to 10 watts when the adverse tilt was increased from 0. 19 cm to 0. 51 cm, during 
post-vibration tests. 

After experiencing a "hard" dryout, the heat pipe did not always recover when the 
load was sharply decreased (up to a factor of ten). As in the qualification testing, re- 
starting the heat pipe after a dryout could always be accomplished by first entering a re- 
verse mode, followed by a controlled' condenser cooldown ( 1° C/MLn), 

Reverse Mode Performance 

Rapid transition from a low load forward mode to a fully blocked reverse mode was 
always readily accomplished by appropriate increases in radiator and reservoir heat inputs. 
All reverse mode tests were run with a melted PCM (>182 R). After a uniform forward 
mode profile was established, a heat load of 40 watts was applied to the main radiator and 
a 4 watt load to the reservoir. The measured response of the evaporator, condenser and 
reservoir temperatures during reversal are shown in Figure 3-21. Both the reservoir 
and condenser temperatures increased in response to the applied loads and after about 2 
minutes, - they both exceeded the evaporator temperature and the evaporator started to 
accumulate fluid, hi 4 minutes, the diode was blocked up to the orifice position. The 
condenser and reservoir were emptied of fluid after 6 minutes, as evidenced by sharply 
increasing temperatures. The evaporator temperature remained constant throughout the 
test while the temperatures of the condenser and reservoir increased by 18° - C and 26° C, 
respectively over an 18 minute' interval. The corresponding temperature rise on the 
blocked side of the orifice was 4° C. 

Transition back to a forward mode was effected by allowing the condenser and 
reservoir to cool and applying a low (2 to, 6 watt) evaporator heat load. 
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Fig. 3-21 HEPP Diode Reverse Mode Response (Post-Vibration) 
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Section 4 


FABRICATION OF A LIQUID TRAP CRYOGENIC THERMAL DIODE HEAT PIPE, 

ENGINEERING MODEL (MOD 9) 

An engineering model of a spiral artery liquid trap diode heat pipe was fabricated at 
Grumman and tested at NASA-ARC. Complete details are presented in Reference 3, but 
the major results are summarized in this section. 

The diode heat pipe consisted of four main parts: an evaporator, a transport section, 
a condenser, and the liquid trap reservoir. The wick was a spiral artery formed by wrapping 
250-mesh stainless steel screen and 0. 04-cm-diam spacers on a' central mandrel. .The 
reservoir had a larger diameter (1.588 cm) than the 0.635 cm OD heat pipe portion and was 
connected to it by a conical transition section and a butt-welded : joint. 1 The internal: core of 
the reservoir consisted of laminated aluminum channels .having a total internal void volume 
of 6.42 cc. There was ho. capillary connection between the reservoir laminates and the 
artery wick. Liquid communication between artery and pipe wall was. achieved with three 
scroll-type webs of 250-mesh screen equally spaced in the condenser and evaporator only. 
Circumferential grooves (63/cm) were used in the evaporator and condenser sections and 
standoffs were provided on the artery for support in the unblocked transport section. 

An overall layout of the liquid trap diode is shown in Figure 4-1 and additional con- 
struction details are given in Table 4-1. . . 

Tests were run in a thermal vacuum chamber with .varying charges of ethane fluid 
(see Reference 3). The optimum charge was found to be 2.67 grams'. The maximum 
measured transport capacity was 12 w-m at 200K with a forward conductance of 6 watts/ 0 C. 
The shutdown energy was estimated as 0.32 watt-hours and the reverse mode thermal 
conductance was calculated as 0.037 W/°C. 

Although the heat pipe demonstrated good transport capacity and reasonable shutoff 
times, it could not be reliably primed either after a thermay dryout or after a cooldown from 
ambient. Reliable startup was only initiated when preceded by a diode reversal followed by 
condenser cooling at a rate less than about 1° C/Min. 
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Table 4-1 Liquid Trap Diode Heat Pipe Design Data 


LENGTH, CM ! 

EVAPORATOR 

10.160 

TRANSPORT 

28.257 

CONDENSER 

30.48 

RESERVOIR 

15.56 

TRANSITION 

1.524 

CHARGE TUBE 

4.445 

INACTIVE SECTIONS 

2.540 

EFFECTIVE 

48.57 

DIAMETER, CM 

PIPE OD 

0.635 

PIPE ID 

0.493 

ARTERY OD 

0.300 

SOLID TUNNEL OD 

0.051 

RESERVOIR OD 

1.588 

CHARGE TUBE OD 

0.318 

OTHER PERTINENT DESIGN INFORMATION INCLUDES. 

PIPE: 

304-1/8 HD STAINLESS STEEL 

SCREENING: 

250-MESH 304 STAINLESS STEEL 

CIRCUMFERENTIAL GROOVES: 

63/CM (160/IN) 

RESERVOIR: 

6061 ALUMINUM LAMINATES 
0.239 .CM THICK WITH 0.127 CM 
DEEP AXIAL MACHINED GROOVES. 

CORE MACHINED TO 1.448 CM OD FOR 
PRESS FIT INTO 304 - 1/8 HD STAINLESS 


STEEL CYLINDRICAL SHELL 

j 0787-038 | 
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Section 5 


THERMAL VACUUM TESTING OF THE ORIGINAL' LIQUID BLOCKAGE (ORIFICE) 
CRYOGENIC THERMAL DIODE, ENGINEERING MODEL (MOD 11) 

( ' K * 

This section presents the results of a thermal vacuum test which was run to check 
previously measured data that had keen obtained using an environmental cold box. The 
details of the earlier, test evaluation are contained in References 2 and 4, The thermal 
vacuum tests confirmed that the performance was similar to the earlier test (November 
1975 and March 1976) and close to theoretical predictions. 

The engineering model liquid blockage diode consists of four sections: an evaporator, 
a transport section, a condenser and an excess fluid reservoir, which is attached to the end 
of the condenser, A blocking orifice is positioned in the transport section about 10. 16 cm 
from the entrance to the evaporator. The liquid reservoir has a larger diameter than the 
rest of the heat pipe (1. 588 cm vs 0. 635 cm) and contains a matrix of laminated aluminum 
channels which retain the excess liquid. There is no capillary communication between the 
reservoir insert and the piral artery heat pipe wick. Design details are summarized in 
Table 5-1 and the overall layout, including thermocouple locations, is shown in Figure 5-1. 

For this, retest, the engineering model cryodiode was configured and instrumented in 
the same manner as in previous tests .(see Reference 4), with one notable exception - the 
charge valve was removed and replaced with an intermediate pinch-off charge tube that was 
joined to the main charge tube by a swagelo c fitting. This minimized the evaporator end- 
effects introduced with a- large mass (the charge valve) and eliminated the need for guard 
heating and cooling. The heat pipe was charged with 4. 8 grams of processed ethane fluid, 
which corresponds to about a 7% overfill at 200 K. 

As before, - three aluminum masses were attached to the pipe: at the evaporator, a 
0. 197' Kg mass to simulate a detector block; a condenser mass of 0. 765 Kg for ease of 
mounting the pipe assembly to an LNg sink; and a reservoir mass of 0.205 Kg, also for 
maintaining the reservoir in contact with' the LN 0 sink. The reservoir and condenser 
blocks were mounted to an LN^ sink which was designed to capture the blocks in a com- 
pression-type arrangement. Rod heaters integral to the compression block were used to 
control the forward mode temperature and also provide an alternate means for diode 
reversal. Strip heaters were attached to the evaporator block for forward mode heat loads 
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Table 5-1 Engineering Model Liquid Blockage Thermal Diode Design Data 


LENGTH, CM 


EVAPORATOR 

10.160 

TRANSPORT-BLOCKED 

10.160 

TRANSPORT-UNBLOCKED 

24.130 

CONDENSER 

21.273 

RESERVOIR 

6.033 

TRANSITION 

1.524 

CHARGE TUBE 

4.445 

INACTIVE SECTIONS 

2.540 

EFFECTIVE 

50.010 


DIAMETER, CM 


PIPE OD 

0.635 

PIPE ID 

0.493 

ARTERY OD 

0.300 

SOLID TUNNEL OD 

0.051 

RESERVOIR OD 

1.588 

CHARGE TUBE OD 

0.318 


OTHER PERTINENT DESIGN INFORMATION INCLUDES: 


PIPE: 304 — 1/B HD STAINLESS STEEL 

SCREENING: 250-rMESH 304 STAINLESS STEEL 

CIRCUMFERENTIAL GROOVES: 63/CM (160/IN) 

RESERVOIR: , 6061 ALUMINUM LAMINATES 

0.239 CM THICK WITH 0.127 CM 
DEEP AXIAL MACHINED GROOVES. 

CORE MACHINED TO I .448 CM OD FOR- 
PRESS FIT INTO 304 - 1/8 HD STAINLESS .' 
' ’STEEL CYLINDRICAL SHELL' 


ORIFICE 

HEIGHT: 

OD: 

LOCATION: 


• .076 CM (.030 IN) 

.483 CM (.190 IN) 

10.16 CM FROM EVAPORATOR ENTRANCE 
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Fig. 5-1 Liquid Blockage Thermal Diode, Engineering Model Overall Layout and Instrumentation 
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and also to the reservoir block as the primary means to initiate diode reversal. All 
thermocouples were copper/ constantan and were welded to the stainless steel pipe envelope. 
They were monitored with a Bristol strip chart recorder. 

The entire assembly was insulated with 25 layers of aluminized mylar super-insula- 
tion (MU) and installed in the thermal lab thermal vacuum chamber. During the tests, the 

-5 

cold wall was not filled and the vacuum was held to at least 10 torr. Thus the cold wall 
was warmer than the operating pipe, which promoted a conservative heat gain estimated at 
less than 0. 15 watt. 

Results: 

Two basic performance points were obtained for comparison: forward mode 
at 1 cm adverse tilt and reverse mode with no tilt. During the forward mode runs the pipe 
temperature was kept near 200 K by suitable adjustments to the condenser block heaters. 
Once the 1 cm tilt was set the evaporator power was incremented and held for 15 minutes, 
until a dryout was achieved. As seen in Figure 5-2, the heat pipe transported 22.4 watts 
, (11. 2 w-m) and dried out when the power was increased to 27 watts. This agrees well with 
.predictions and also indicates that performance had not degraded from the earlier tests. 

•The forward mode conductance as calculated from the measured thermal profiles was 
3.6 w/°C (see Fig. 5-3), 

During reverse mode testing, the reservoir heater burned out and it was necessary 
to initiate reversal by overpowering the condenser block heater. The evaporator load was 
zero. As observed in Figure 5-4, a transient plot of average temperatures, the condenser 
temperature increased at a rate of 1.4° C/min and the temperature rise of the reservoir 
preceded that of the condenser due to its smaller mass. The evaporator temperature, 
which .started at -79° C, held virtually constant and reached -76° C at the end of 30 minutes 
when the condenser temperature registered -51° C. Full evaporator blockage was achieved 
about 10 minutes after starting. 
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Fig. 5-2 Cryodio'de Test Points 
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Fig. 5-3 Test Temperature Profile for Engineering Model Liquid Blockage Thermal Diode 
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Fig. 5-4 



During Reverse Mode 







Section 6 


: AN INVESTIGATION OF THE START-UP DYNAMICS OF SPIRAL. 

ARTERY CRYOGENIC THERMAL DIODE HEAT PIPES (MOD 12} 

The .previously reported test results for the cryogenic thermal diodes were good ex- 
cept for one persistent problem. The only way they could be reliably started from a dried 
but condition was to heat the condenser, which flooded the evaporator section with ethane 
fluid, followed by a gradual transition to a forward mode thermal gradient (evaporator 
temperature above condenser temperature). This start-up technique proved necessary for 
both the liquid trap and liquid blockage- designs. Attempts to start the pipes directly froni 
a room temperature cooldown or from an adiabatic condition after a hard dry-out usually ' 
failed, although a few successes were witnessed, usually after an initial cooldown; 

Presuming a similar characteristic in zero g, this means that if a cryogenic thermal 
diode is used as a piece of flight hardware there would have to be operational constraints ' 

s ( 

imposed which would involve preconditioning the heat pipe by flooding the evaporator. 
Within the present state-of-the-art, as far as a flight application is concerned, this means 
that only one preconditioning reversal mode is necessary as long as the imposed thermal 
loads are below the dry-out limit. In actuality, this might be acceptable for most cases 
since design safety margins of 2 to 4 are typically used. However, in some cases even the 
one preconditioning sequence might not be acceptable. Also, there is no reason to believe 
that this start-up problem is peculiar to diode-artery pipes, it probably exists in any cryo- 
genic artery-type pipes made with the same materials and fluids. Since there are many 
useful spacecraft and sensor applications for cryogenic heat pipes and thermal diodes, it is 
important that they can be used without imposing any operational constraints upon the 
mission. ' . 

The specific reason for the start-up problem is not clear, but there are two prime 
candidates: ' surface wetting and artery priming. A surface wettingproblem could exist 
between the 304 stainless steel surfaces and the ethane fluid, hi a dried out condition, the 
existing surface chemistry (possible as a result of contaminants) could cause a non-wetting 
situation which can only be overcome by forcibly flooding and hence conditioning the sur- 
faces. Once the surfaces have been wet, they tend to remain wet until a dryout condition 
is experienced. The "harder" the dry-out, the more surface area becomes dried out and 
non-wetting. Thus, on a subsequent start-up attempt the active evaporator surface area is 
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severely restricted which results in much lower throughput capability. There is strong 
evidence to support this wetting theory, since once a heat pipe has been Started it can be 
repeatedly cycled from no load to high load as long as a dry-out situation is not encountered. 
]h fact, it appears that, once started, it can remain in an adiabatic. state indefinitely before 
a load is applied. 

The other theory explaining the difficult start-up is that an artery priming problem 
exists. This would most probably be caused by the entrapment and subsequent growth of 
vapor bubbles within the arterial voids, thereby interrupting the liquid flow. Consider a 
cool-down from room temperature. • Initially, the entire system (vapor and liquid) is iso- 
thermal and in an adiabatic condition. The screen artery is partially saturated by virtue 
of the. screen-wieking capability. During the cool-down, heat is being pulled out of the 
condenser and the vapor pressure can drop below the saturation pressure of the liquid in 
the artery. -The warmer liquid within the artery then tends to flash, which results in the . 
formation of internal' vapor bubbles. At this time the evaporator is dried out because it 
cannot draw enough fluid from the vapor-blocked artery. Gradually as the artery liquid 
cools, the saturation pressure drops and liquid capillary forces tend to collapse the 
arterial vapor bubbles, hi theory, a partially primed artery can be fully primed by per- 
mitting continual fluid circulation at low heat loads and the introduction of cooler liquid 
from the condenser to the other parts of the artery. After a time, the cool liquid advances 
to the evaporator where the last remants of vapor bubbles collapse under the capillary 
forces. This sequence can be accomplished by operating the deprimed heat pipe at a load 
which is less than the maximum capacity of the partially primed artery and keeping it 
operating long enough to collapse all of the yapor bubbles. 

Priming problems that come after a hard dryout at operating temperature can be 
described as follows. Once the capillary pumping limit is exceeded, the amount of fluid 
vaporized decreases which subsequently causes the vapor pressure to drop. This results 
in a pressure imbalance across the liquid in the artery which causes liquid flashing and 
vapor bubble formation within the artery. The displaced liquid ends up in the heat pipe as 
free liquid and the deprimed artery has a much reduced load, capacity. Once again, in theory, 
the artery can be primed if operated long enough at a low load within the capability of the 
deprimed wick; As before, this is accomplished by the cool condenser liquid causing pro- 
gressive collapsing of the vapor bubbles as it advances toward the evaporator. 
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For both the depriming sequences cited, the artery can be forcibly reprimed by 
heating the condenser. This raises the vapor pressure which then forces free liquid to 
enter the artery and displace or collapse the entrapped vapor bubbles. 

The actual reason for the stari-up difficulties encountered with the ethane artery 
cryodiode heat pipe. could well be a combination of both wetting and priming problems. 
However, it was felt that the greatest insight could be'obtained by separate investigations. 
Surface wetting characteristics were studied by experimental means, and an analytical 
investigation of the priming phenomenon consisted of predicting detailed temperature dis- 
tributions and their influence on bubble growth. 

6.1- ANALYTIC STUDIE S OF PRIMING 

Introduction 

The 1/4" (6.35 mm) OD cryogenic diode with spiral artery wick has been made. in 
several versions and tested with a variety of working fluids. Many of the diodes experienced 
difficulty in priming and/or required special priming procedures, such as diode reversal, 
prior to successful forward mode operation. To provide a better understanding of the artery 
priming problem, an analysis was performed which takes into account the detailed internal 
structure of the artery. The gap regions of the artery are assumed to be initially filled with 
vapor. The artery is divided longitudinally into sections, and calculations are performed to 
determine if each section, starting at the condenser end of the pipe, can fill. If all sections 
fill, the artery is said to be primed. 

The discussion of the analysis covers the following areas: 

1) Heat pipe geometry being studied 

2) , Priming criterion 

3) Nodal model 

4) Calculation of thermal couplings used in nodal model 

5) Computer program 

6) Analytic results 

7) Conclusions and recommendations. 

While only one pipe configuration was evaluated in this study, the method of analysis is 
easily adapted to other geometries and working fluids. 
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Diode Heat Pipe Geometry 

The diode geometry modeled in this analysis is shown schematically in Figure 6-1. 

The pipe has a 1/4 inch (6,35 mm) OD and, exclusive of reservoirs, is 27 inch (0.68 m) 
long. For purposes of this study the reservoir is neglected. A liquid blockage reservoir 
would be the first region to fill on cooldown, and would be inactive during the priming 
process. A liquid trap reservoir would-be dry during initial cooldown, and also be inactive 
during the priming process. 

The pipe has a 4 inch (10. 16 cm) evaporator, 15 inch (38. 1 cm) adiabatic section, and 
8 inch (20.32 cm) condenser. The evaporator uses electrical strip heaters for heat source 
and the condenser is clamped to a liquid nitrogen cooled block. Evaporator and condenser 
are wrapped in multilayer insulation (MLI) and the complete assembly is assumed to operate 
in a vacuum environment. 

Priming Criterion 

The priming criterion used in this analysis is that a vapor bubble in the gap region 
of the artery must collapse when subject to a temperature distribution within the pipe which 
is characteristic of steady state heat transfer conditions. 

In equation form, this requires 

P.v, g < PI, g + A Pc, g 

where Pv, g is the saturated vapor pressure characteristic of a bubble within the artery, - 
PI, g is the liquid pressure adjacent to the bubble at the point on the bubble, closest to the 
condenser and highest in elevation in a 1-g field, and APc,g is the pressure difference 
across the meniscus (surface tension divided by half the gap width) at the condenser end of 
the vapor bubble in the artery. The vapor bubble is assumed to be in thermal equilibrium 
with the liquid filled gap and screen regions which surround it. This implies some transport 
of heat (as vapor flow) within the bubble. There will also be an associated (extremely small) 
temperature variation over the surface of the bubble, which is neglected in this analysis. 

The bubble is treated as isothermal. 

At the onset of this priming study, some thought was given to performing a transient 
calculation, modeling the pipe in detail as cooldown proceeded, starting from an initial iso- 
thermal room temperature condition. There were several difficulties with this approach: 
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1) Some fluids (e.g. methane) would initially be above critical temperature, presenting 
modeling problems for beat transfer within the pipe. Below critical temperature, 
there would be a problem in handling the filling of the porous regions with liquid 
eonsistant with the structure of a nodal thermal analysis model, and computational 

. stability problems associated with the vapor nodes (small thermal capacitance and 
large thermal couplings). 

2) Transient cooldown curves can vary considerably, so that a number of curves 
would have to be evaluated, 

3) Programming and computation time would be considerable compared with a steady 
state approach. 

Difficulties aside, however, it was recognized that in most cases where priming- 
problems were encountered, the experimentalist would let the pipe remain (soak) at cold 
conditions for several hours. Hence, the ability of the pipe to prime under cold steady 
state condition is, in a sense, a more important aspect to assess analytically. ■ A pipe 
which failed to prime under these conditions would not be likely to prime during a mono- 
tonic cooldown. 

Nodal Model 

The nodal model is shown schematically in Figure 6-2. The pipe is divided longitudin- 
ally into seven sections (2 in the evaporator,. 3. in the adiabatic section, and 2 in the con- 
denser). Multilayer Insulation (MLI) surrounds the evaporator and adiabatic sections. 

The condenser is coupled directly to a cooling block. In the radial direction, the MLI 
(which is subject to large temperature differences) is divided into 2 nodes, and the pipe 
itself into 4 nodes - the pipe wall, the vapor space, the porous wicldng (artery screening 
and webs) and the gap region within the artery. A description of the nodes follows: 

1) Nodes 1-5 - outer layers of MLI (nodes 1 and 2 cover the evaporator and nodes 
3-5 cover the adiabatic section) 

2) Nodes 6-10 - inner layers of MLI 

3) Nodes 11-17 - pipe wall 

4) Nodes 18-24 - pipe vapor space 

5) Nodes 25-31 - porous wieldng consisting of three cylindrical webs and the screening 
used to make the artery 
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Fig. 6-1 Diode Heat Pipe Geometry 
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Fig. 6-2 Nodal Model for Priming Analysis 
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6) Node 32 - the vapor filled gap region within the artery, initially, node 32 extends 
the full length of the artery and is thermally coupled to all the surrounding porous 
nodes. Priining is then tested for a meniscus across the gap at the condenser end 
(node 31). If priming can succeed, the gap region within node 31 is assumed filled 
with liquid, and the coupling between nodes 32 and 31 is set equal to zero. This is 
equivalent to having the vapor region of the gap extending from nodes 25 to 30. 

A new steady state temperature distribution is then obtained, and priming tested 
for filling node 30. The process is repeated until either priming fails or the last 
node (node 25)is filled 

7) Node 101 - A boundary node representing the inside wall of the vacuum chamber, 
assumed kept at a constant temp (90° F) (27° C) in this analysis. 

8) Node 102 - A boundary node representing the liquid nitrogen coolant in the con- 
denser block, assumed to be at -136° F (-93° C) in this analysis. 

A brief description of the couplings is given below. A more detailed description 
showing the quantitative evaluation of each coupling is presented in the next section. 

t 

A total of 22 radiation couplings are provided: 

1) Vacuum chamber to outer MLI (5 couplings) 

2) Outer MLI to inner MLI (5 couplings) 

3) Inner MLI to pipe wall (5 couplings) 

r * 

. .4) Pipe wall to porous material (7 couplings) 

It may be noted that no radiant couplings have been provided in the longitudinal direction. 
Surrounding the MLI, the vacuum chamber is considered a relatively large high, emittance 
enclosure. Within the MLI, longitudinal radiant heat transport is included in the couplings 
for conduction parallel to the laminations.' Within the heat pipe, radiant flux is generally 
small because of the low temperatures and the view factors for longitudinal transport are 
fairly small. 

A total of 43 ordinary conduction couplings are provided: 

1) Parallel conductance within MLI (8 couplings) 

2) Heat pipe wall (6 couplings) 

3) Condenser wall to LN^ (2 couplings) 

4) Evaporator wall to vapor (2 couplings) 


6-8 



5) Adiabatic wall to vapor (3 couplings) 

6) Condenser wall to vapor * <2 couplings) 

7) Vapor to porous wi eking (7 couplings) 

. 8) Porous wicking to gap region (7 couplings) 

9) Between vapor space nodes (6 couplings) 

A total of ten fluid conductors (one-way couplings) are provided; 

1) Condenser wall to porous wicking (2 conductors) 

2) Within porous wicking (6 conductors) 

. 3) Porous wicking to evaporator wall (2 conductors) 

It may be noted that fluid conductors are provided for the liquid flow path, whereas 
ordinary conductors are provided-for the vapor flow path. The use-of ordinary conductors - 
for the- vapor flow is related to the computational procedure- used to calculate the vapor 
flow. For each vapor node, the vapor is assumed to be- at saturation pressure, and the 
vapor flow between nodes is a consequence of the node to node pressure differences, - 

Thermal- Coupling Calculations (Engineering Units' used for convenience) 

Nodal Geometry 

A schematic of the nodal geometry, showing lengths and radii, is presented in Figure 
6-3. Nodes are taken as centered within the boundaries. 

Radiation Couplings (22) 

1) MLI Blanket 

The MLI blanket is assumed to consist of a total of 25 layers of crinkled, single 
aluminized mylar. This is divided into inner and outer blankets, each having an effective 
emittance of 0,015, Then 


Ri,j ^gff 5 “ ^ e ££ 27T r i,j li 

For example 

Rl,6 = 0.015 x 27 T 

The coupling between the outer MLI nodes and the vacuum chamber, R^ 1 increased 
to 6,545 x 10~ 4 to reflect only half as many layers, and this value was also used between the 
inner MLI node and pipe outer wall, R g 11 . The corresponding blanket effective emittance 
values, allowing for the change in mean radius, are 0.0166 and 0.0436, respectively. Some 
runs were made with both of these values taken as 0.024. Results were similar. 
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2) Pipe Inner Wall to Artery Surface 

The coupling between the pipe inner wall and the artery outer surface is based on an 
effective emittance for both surfaces of 1. 0 and an area based on a mean radius of 0. 153 
inch less the thickness of the three webs (3 x . 041 inch) which support the artery. 

Table 6-1- is a list of the radiator couplings, showing nodes coupled and the coupling 
values. 


Conduction Couplings - (Ordinary) 


1) Longitudinal Conductance in MLI (Parallel to the Laminations)' 

Thermal conductance within the MLI parallel to the laminations is a complicated 
radiation and conduction phenomena which does not lend itself to treatment as a pure radia- 
tion or pure conduction coupling. Experimental data is available (Ref. 5) for crinkled 
single side aluminized mylar expressed in terms of effective conductivity, K ^ f , at layer 
densities of 1850 and 3230 layers/meter. The data is shown plotted in Figure 6-4, from 
which it can be seen that is a strong function of temperature (taken here as a fourth 
power mean temperature). 


We have 25 layers of MLI in a 3/8 inch thick annulus surrounding the pipe, correspond- 
ing to a layer density of 2625 layers/meter. An estimated curve for this layer density is 
also shown in Figure 6-4. The curve is a straight line on the semi-log plot with the form 


where T. . = 

*-9 J 


K „ = a f b 
eff 

4 4 “ 

T. +T. 
i J 




1/4 


a = 1.329 (10 ) w/(m-k) 

b = 1.972 


the conductance coupling is then 


K. . = K w 
i,j eff 


x 


^»3 


where A is the cross-sectional area for parallel conductance 
x 


x 


7 2 2 , 

- W <r D - rj ) 
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Table 6-1 Radiation Couplings 


NO. 

' NODES COUPLED 

B, (FT 2 ) 


1 

1,6 

4.091 x 10" 4 


2 

2,7 

4.09 x 10 4 


3 

3,8 

10.23 x 1 0’ 4 


4 

4,9 

10.23 x 10 4 


5 

5,10 

10.23 x 10 4 


6 

101,1 

6.545 x 10" 4 


7 

101,2 

6.545 x 10 4 


8 

101,3 

16.36 x ID" 4 


9 

, 101,4 

16.36 x 1 0' 4 

1 

10 

101,5 

16.36 x ID" 4 


11 

6, 11 

6.545 x 10 4 


12 

7, 12 

6.545 x 1 0 -4 


13 

8, 13 

16.36 x 1 0' 4 


14 

9,14 - 

16.36 x 1 0' 4 


15 

10, 15 

16.36 x 1 0' 4 


16 

11,25 

4,968 x 1 0' 3 


17 

12,26 

4.968 x 1 0' 3 


18 

13,27 

12.42 x 10' 3 


19 

14, 28 

12.42 x 10~ 3 


20 

15, 29 

12.42 x 10‘ 3 


21 

16, 30 

9.935 x 10~ 3 


22 

17,31 

9535 x 1 0’ 3 
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Fig. 6-3 Schematic for Nodal Geometry 


6-11 






























Fig. 6-4 for Crinkled Aluminized Mylar 
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and 1. . is the distance between node centers. The A values for the outer and inner 
^5 J X 

blankets are 3. 32xl0~ 3 ft 2 and 1. 79xl0“ 3 ft 2 respectively. 

Because of the strong temperature dependence, the K. . values must be computed 

l >3 

based on nodal temperatures, making the solution interactive. The initial values loaded in 
are based on a mean temperature of 200° K (360°R). At 200° K, = . 04583 w/m- °K 
(.0265 BtuAr ft °R) 

Then, in English units. 


K. . = 0.0265 




Values for K. . are given in Table 6-2 for 360°R, for the 8 longitudinal MLI couplings. 
2) Pipe Wall Conduction Couplings 


The pipe wall is stainless steel with an assumed thermal conductivity, k - 10 Btu/ 
(hr-ft - °F). The couplings' are. then , , 


..K. . = KA = K 

,x „ 


7T (d o 2 -d. 2 ) 


L ijj 


L. . -4 

1,1 


where d = 0.250 inch and d. = 0.194 inch, and L. . is the longitudinal distance between 
node centers. 


Resulting K. . values are given in Table 6-3,- 
3) Condenser Cooling Block Conduction Couplings 

A schematic of the condenser cooling block geometry is shown in Figure 6-5. The 
total thermal resistance is made of of 5 parts - resistance between the node center and the 
pipe -surface, contact resistance between pipe and V-Block, resistance across V-block, 
contact conductance between V-Block and LN^ cooled block, and resistance across the LN^ 
cooled block. Each node has a 4" longitudinal length. Then 



JL 

L 


t/2 

k 7 r S 

ss 




- 1 , 1/2 • 1/2 

7TD h . S K ^ yh 
o o v al - o 



x 


6-13 



Table 6-2 MLl Parallel Conductance Couplings (T = 360°R) 



Ky I at 360°R) 
BTU/HR °F 


5.285 x 10 
3.020 x 10"* 
2.1 14 x 10^ 
2.1 14 x 1.0 J 
2.846 x 10" 
1.626 x 10^ 
1.138 x 10 J 
1.138 x 10^ 
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where L = node longitudinal length - 4 inch 

' t = pipe wall' thickness - .028 inch 

D"= * pipe node to wall mean dia - 0.236 inch 

k = conductivity of stainless steel - 10 BAr ft ° 
ss 

D q = pipe OD - 0.250 inch 

2 

h Q = contact conductance = 1000 BAr ft 0 F 

S' = shape factor, V Block = 0.9601 

k^ = conductivity of aluminum = 100 BAr ft 0 F 

y = height of V-Bloek contact = 1. 375 inch 

x = heat transfer length in LN^ block = 1 inch - 

Then, upon substitution of the appropriate values, 

1 = 1 

K. .• 11.0 

. . 

or 


■ ■ K i,j K 102,16 

These are conductors 15 and 16. 


K 102,M = 11 - 0B/hr ” :P 


4) Pipe Wall to Vapor Coupling’s 


The measured film coefficients for evaporation and condensation for ethane at 200 K 

2 

are 1620 and 2890 BtuAr ft °F, respectively. The conductance values are then 


K i,i “ h7r Vm 


e.g., 

/ 0.194 \ 2 

K 11,1S - 1620 x * (“IT ) 5 “ 13 ' 7 BtaAa • 

For the adiabatic section, the wall will in general be hotter than the .vapor, and dry, 

since there is no wicking between wall and artery in this region. For laminar vapor flow, 

2 

the film coefficient is estimated to be 11 BtuAr ft °F. 


The resulting wall to vapor couplings are given in Table 6-4. 
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5) Vapor Space to Artery Couplings 

Evaporation or condensation may take place at the artery surface, with film co- 
efficients representative of the thermal resistance of the liquid layer between the artery 
node and the artery surface. The conductance is then 



77-15 f 1* f 

art 1,3 


The thickness of the liquid layer is estimated as one layer of screening (since the outer 
sheath of the artery is 2 layers thick). Then, e.g. , 

7T(0, 112) 2 

K 18,25 ~ 0,0797 0.0026 12 = 1,798 BtuAr °F 


The set of resulting values is given in Table 6-5. 

6) Artery Gap to Screen Couplings 


Within the artery, the coupling between the vapor filled gap regions (node 32) and 
the spiral screening is based on the thermal resistance of one layer of screening. The 
area is based on a gap width inside the artery of 0. 318 inch (assuming only the outer surface 
of the gap is effective for heat transfer), and the length of each porous screening node. 



w 1. 

l 

t 


e.g.» - 

0.318 (2) ‘ 

K 9 , = 0.0797 x = 1.635 B/hr °F 

9 ‘ (.0026) (12) 

The set of resulting values is given in Table 6-6. 

7) Vapor Flow Couplings 

For any vapor node, the vapor pressure is set equal to the saturation pressure 
corresponding to the temperature of the node. The flow rate of vapor between any two 
adjacent nodes, assuming laminar flow, is then 
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Table 6-3 Pipe Wall Conduction Couplings 


NO. 

NODES COUPLED 

Kjj BTU/HR °F 

9 

11,12 , ' - 

8.137 x 10" 3 

10 

12,13 

4.650 x 10‘ 3 

11 

13,14 

3.255 x 10" 3 

12 

14,15 

3.255x 1 O' 3 ” 

13 

15,16 

3.616 x 10' 3 

14 

16,17 

4.068 x 10' 3 
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Table 6-4 Wall to Vapor Couplings 


NO. 

NODES COUPLED 

Ky BTU/HR °F 

17. ... .. 

; . . 11, 18 _ ' , 

13.7 

18 

12,19 

13.7 

19 

13, 20 

0.236 

20 

14,21 

’ 0.236 

21 

15,22 - - * 

. ’ • . 0.236 

22 

16,23 

48.9 

23 

17,24 

48.9 
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Table 6-5 Vapor to Artery Conductance Couplings 


NO. 

NODES COUPLED 

Kjj BTU/HR °F 




24 

18,25 

• - 1 .798 

25 

- 19,26 

1.798 

26 

20, 27 

4.494 

27 

*' '21.28- ‘ 

4.494 

28 

22, 29 

4.494 

29 

23,30 

3.595 

30 

24,31 

3.595 
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A D 2 / P G 

V V / VC 


m. 


l * J 2 (f N r )\ fx. 


'P - P 
1 


where A = vapor flow area 

v 

D v = hydraulic diameter of vapor space 
f = friction factor (shear stress/dynamic head) 

N,, = Reynolds Number 

±t 

P v = Mean vapor density for nodes i, j 
- Mean vapor viscosity for nodes i,j 

= gravitational constant 

P P 

v i’ Vj = vapor pressure for nodes i, j 

L. = distance between node centers for nodes i, j • 

i. - 

3 

As written, m. . is positive for flow from i to j and negative for the reverse. The 
i» J 

heat flux associated with this mass flux is 

Q. . = h„ m. - 

t> 3 f g 1,3 

where h. = mean latent heat for nodes i,j 
g 

In terms of conductance, K. . 

i,3 

Q. . = K. . (T. - T ) 

1,3 1,3 i j 

The pressure and temperature differences are generally very small. To avoid numerical 
accuracy problems, the approximate form of the dausius-Clapeyron relation is used: 


P - P 
V i V j 


'd P 


v 


P h» 
v fg 


dT /i,j 

where J = mechanical equivalent of heat (778 ft-lbs/Btu) 


T. ' - T. 
i 3 


T. . 
1,3 


and T.. = (T: +T.)/2 
i] i y 
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Also defining 


k. . 
M 


A D " 
v v 


S„ 


2(EN r ) L : 


i,j 


we can-write: 


K. . 


= k, 


Uj 



Values for k. .■ are input as initial values for the K. . couplings. The CL Cl subroutine saves 
i>3 • L3 

all the initial conductor coupling Values, and recomputes new K. . values for the latest set of 

. • * LJ 

temperatures in the nodal solution. 

In evaluating the k. . , the vapor space area is the annullar area between wall and artery, 
i» J 

less the web cross-sectional area, the hydraulic diameter is taken as twice the annulus area, 
and the f3SL> product for laminar flow is taken as 23; 8. Then evaluating, for example, 

K 


18,19* 


(. 01575/144) (. 082/12) 2 (4. 18xl0 8 ) 


"18, 19 


2 (23.8) 


( 2 / 12 ) 


= 0.2691 


lb - ft 
m 

lb. - hr 2 
1 


Values for the k. . are given in Table 6-7. 
i» J 

Liquid Flow Conduction Couplings - ■ ■ 

The liquid-couplings are one-way conductors, identified in the nodal model by inserting 
a minus sign in front of the upstream node number. The heat conducted may be written 


whence 


Where 


Q . . = K . . (T. - T ) = m . . 
-i,J -1,3 i o' -1,3 

K . . =■ m . . C_ . 

-i,j . -1,3 Pj- - 

C-p = liquid specific heat 
*1' 




The liquid mass flow is determined from the vapor mass flow, since the liquid inventory at 
each node is constant for steady state conditions. For example: 


m 


-25, 11 


m 


11,18 


That is, the liquid flow from the artery node 25 to evaporator node 18 must equal the 
evaporation at evaporator node 18. 
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Table 6-6 Artery Gap to Screen Couplings 


NO. 

NODES COUPLED 

Ky BTU/HR °F 

31 

25, 32 

1 .635 

32 

26, 32 • 

, 1.635 

33 

. 27,32 

4.07 

34 

28, 32 

4.07 ' 

35 

29,32 

• 4.07 - ‘ ' 

36 

30,32 

3.27 

37 ' 

31,32 ‘ 

3.27 ’ • 

0787-045 




Table 6-7 Vapor Passage Couplings; K, ; 

'/J 


.* , CONDUCTOR NO. 

NODES COUPLED 

K. . (LBM FT 4 /LB HR 2 ) 

I J 

' ‘ ’ '38' 

39 

40 

41 

42 

43 

0787-046 

18.19 

19.20 

20. 21 
21,22 
22,23 
23 24 

0.2691 

0.1538 

0.1076 

0.1076 

0.1196 

0.1345 
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Then since 



K (52) K -26 , 12 ~ C 2 K 12,19 ^ T 12 " T 19* 

For flow along the artery from node 26 to node 25, conservation of mass at node 25 gives: 

K < 51 > = K -26, 25 = K -25, 11 + °2 [ K ^ < T 25 * T 18> + K < 31 > < T 25 ' T 32> 

where the extra terms account for evaporative flow to the vapor space and gap. Similar 
expressions conserving mass at nodes 26, 27, 28, and 29 (in sequence) are: 

K < 60 > = K -27, 26. “ K < 5l > + K < 52 > + C 2 [ K < 25 > < T 26 ’ T 1S> + K < 32 > < T 25 ' T 32>] 

K < 49 > = K -28,27 ° K < 50 > + C 2 [ K < 26 > < T 27 * T 20> + K < 33 > <*28 " T 32>] 

K < 48 > = K- 29,28 ’ K < 49 > + C 2 [ K (27) < T 28 " T 21> + K (34 > (T 27 ' T 32>] 

K («) - K -30,29 - K < 48 > + C 2 [ K < 28 > < T 29 - T 22> + K < 35 > <*28 ' T S2>] 


The liquid flow from the condenser wall to the artery is assumed to he the same for 
each condenser' node (since the condenser nodes could exchange liquid between themselves), 


giving 

K < 45 > * K -16,30 * — [ K < 22 > < T 23 - T 16> + K < 23 > < T 24 ' T 17>] 
and K (44) = k _ 17j31 = K < 45 ) 


Then a mass balance at node 30 gives 

K <«> - K-31, 30 " K < 47 > ' K < 45 > + °2 [ K < 29 > < T 30 - T 23> + K < 36 > ^SO ' T 32>] 

These couplings are recalculated with each iteration in the nodal model solution. 

Initial values of 1. 0 are input for each liquid flow coupling for data handling purposes, but 
these are replaced by computed values during the first iteration. 
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Computer Program 


The computer program used in this analysis is a large general purpose nodal program 
called T-15 - Transient Thermal Analyzer. T-15 contains a number of subroutines, in- 
cluding STEADY, which iteratively determines temperature distributions at steady state. 
T-15 also provides subroutine CLCI called during each iteration and subroutine CLC2, 
called after determination of a steady state solution. These subroutines can be used to 
update the nodal model or perform any auxiliary operations desired by the user. 

For this artery priming study, the CLCI subroutine was used to perform the following 
operations: 

1) Store the initial data set of conductance coupling values, for possible use later in 
the calculation 

2) Modify couplings as required for parametric studies (vapor flow resistance, MLI 
conductance parallel to the lamination, and MLI effective emittance) 

3) Calculate latest values for temperature dependent couplings (MLI parallel con- 
ductances, MLI normal emittance, vapor conductance, fluid conductors) during 
the iteration for the steady state temperature distribution 

4) Calculate the priming head at the condenser end of the artery, with a statement 
"Artery primes in 1-G" printed out. This is not the real test of priming (which 
is made in subroutine CLC2), but it does provide evidence of each pass through 
CLCI during the calculation. 

The operations performed by the CLC2 subroutine are: 

1) Verification that priming can start at the condenser end of the artery 

2) Determination of artery permeability for liquid pressure drop calculations 

3) Modification of the nodal model to reflect filling of the artery node at the condenser 
end of the pipe (node 31) with liquid. This involves setting the coupling between 
this node and the gap (coupling 37) to zero, and reducing the coupling to the vapor 
space (coupling 30) to reflect the increased liquid thickness. 

4) Recalculation of the steady state solution and printing the result 

5) Calculation of the priming head in 1-g, allowing for the liquid pressure drop in the 
liquid filled artery node, and printing whether or not priming of node 30 can take 
place with node 31 filled 
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6) If priming can take place, the couplings between node 30 and the gap and vapor 
space are modified as in step (3) above to reflect node 30 filled with liquid, and 
the calculations of steps (4) and (5) are repeated. This process is similarly re- 
peated in sequence for priming of nodes 29 down to 25. If node 25 can fill with 
liquid in 1-g, the entire artery is primed 

7) If priming in 1-g cannot occur at any of the nodes from 31 to 25, the gravity head 
(liquid density times the gap height above the bottom of the pipe) is removed to see 
if priming can occur in o-g. 

8) If priming can occur in o-g, the node is assumed to fill with liquid, the couplings 
are adjusted as in step (3) above, a statement is printed out that priming in o-g 
can occur, and steps 4-7 are repeated. If node 25 can fill with liquid in o-g, the 
entire artery is primed 

9) If priming cannot occur in 1-g, following the o-g priming caluclations, the evaporator 
heat flux is reduced (by an amount specified, in watts, as dummy conductor 127) and 
the entire sequence of steps (1) thru (8) is repeated. 

Subroutines CLC1 and CLC2 are .included in Appendix B. 

Analytic Results 

The computer nodal model previously described was mm for off-design variations of 
three parameters; wick permeability C(124), MLI effective emittance C(123), and MLI 
parallel conductivity C(126). Wick permeability was found to have a direct influence on the 
ability of the artery to prime itself when under load. As seen in Figure 6-6, the maximum 
evaporator heat load under which the artery will still prime (i. e. , maximum priming load) 
changes nearly linearly with wick permeability. The nominal artery design (normalized 
permeability = 1. 0) has a maximum priming load of 5 watts, which decreases to under 3. 5 
watts at 75% of the design value. This type of performance is consistent with observed test 
results. Figure 6-6 indicates that higher priming loads can be obtained by suitable in- 
creases in wick permeability. For instance,, doubling the permeability would increase the 
maximum priming load to about 10 watts. Another point of interest is that the o-g priming 
load is only l/2 watt higher than the 1-g value. This means that any artery priming pro- 
blems which exist in 1-g will not vanish in o-g. 

The effect of degradation in MLI performance is to cause the wall temperature in the 
adiabatic (transport) section to become warmer than the evaporator wall temperature. 

This can be explained by the fact that the transport section is more poorly coupled to the 



MAXIMUM PRIMING LOAD, W 





condenser section and is therefore more sensitive to changes in MLI conductance which 
increase the coupling to the surrounding higher temperature ambient. Typical results are 
s umm arized below for degraded values of both effective emittance ( € and parallel 
conductivity ( € using the nominal wick permeability. For these cases, the condenser 
temperature sink was -93° C. 


Normalized MLI Parameter 

(Tevaporator - Tadiabatic)' 

c eff 

*11 


1.0 

1.0 

- .23 

2.0 

1.0 

- .148 

10.0 

1.0 

‘+ . 096 

10.0 

10.0 

+ .423 


As seen, with nominal MLI parameters the evaporator is colder than the adiabatic section 
(At is negative) as would be expected. However, as the MLI degrades the At becomes 
more positive, which means the adiabatic section is warmer than the evaporator. This 
curious occurrence has actually been observed during some thermal vacuum tests, which 
correctly infers that these test articles were poorly insulated. 

Conclusions and Recommendations 

To permit higher maximum priming loads, the wick permeability must be increased. 
This can be accomplished by using composite artery wick designs with coarse interior mesh 
and large flow channels surrounded by the finer screen mesh needed to create the necessary 
capillary pressure rise. Also, it is recommended that permeability measurements be per- 
formed on each artery in addition to the normal capillary rise test which is used to check 
artery integrity. 

6.2 EXPERIMENTAL INVESTIGATION OF SURFACE WETTING 

Laboratory tests have been conducted to determine the wettability of ethane with 
various stainless steel wick samples. This data is considered necessary to help understand 
and explain priming and transport behavior recently obtained with the subject diode heat pipe. 

Ethane wettability was determined by measuring its wicking height at cryogenic 
temperatures and comparing it with other fluids such as acetone and Freon-113 whose 
wetting characteristics are well known. The wicking height to which a particular fluid will 
rise with a square mesh screen wick is given by: 
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h 


2 cr cos 9 

r p 


( 1 ) 


where h - measured wi eking height 
O* = surface tension 
9 = contact angle 

r = pore radius of screen opening, taken as one-half of opening 
p = liquid density 

A fluid which has perfect wettability will have a zero contact angle, and therefore, 
its ideal wicking height, h. is . 

2 <T 

h. =— — 
x r p 

Thus, a measure of wettability of a fluid with a particular wick is its actual wicking height 

versus its ideal wicking height, or 

i m _ 2cr cos0 /r P * 

wettability = h/h. =?- — — ~ cos 0 

2 cr/xp ■ 

• A test rig was designed to allow visual observation of the wicking height (h) which, 
along with known fluid and wick properties ( P, r) was used to calculate wettability or 
cos 6 from equation (1). The wettability test apparatus is shown in Figure 6-7, It consists 
of a welded aluminum housing with three flanged ends that provide a leak tight enclosure 
using O-ring seals. The top flange has a see through glass port for viewing into the chamber. 
The test sample is held in position directly beneath the view port. One side flange has pro- 
vision for evacuating the enclosure, introducing charge, and monitoring internal pressure, 

A cooling coil is wrapped around the chamber to enable it to reach cryogenic temperatures 
when testing with ethane. Thermocouples on the chamber wall monitor system temperature. 

To determine the wicking height, a sample of screen wick about 1 inch (2. 54 cm) by 
4 . -inches (10. 16 cm) is inserted into the chamber and held in place using a small spring 
clamp. The long side of the sample is notched every 1/2 inch (1, 27 cm) to permit estima- 
tion of the wicking height attained. The chamber is then evacuated and cooled to the desired 
temperature. Next, the test fluid is introduced from a specially prepared charge bottle. 

The apparatus is tilted horizontally until the entire sample is wet with fluid and then set 
in the vertical position and allowed to stand until the wetted height reaches equilibrium (see 
Figure 6-8). The bottom of the wick is always in contact with fluid which will rise to a 
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a) TILTHIG HORIZONTAL TO 
WET ENTIRE SAMPLE 


Q7B7-034 



Fig. 6-8 Test Procedure 
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particular height on the wick according to its wettability. By shining a narrow beam of 
light through the viewing port, the darker wetted height can be discerned from the lighter 
unwetted region and its height gaged by the 1/2 inch (1. 27 cm) marks along the sample. 

A description of the wick sample used is given in Table 6-8. Applicable fluid 
properties are listed in Table 6-9 for the three fluids tested. Freon-113, acetone and 
ethane. 

Test restilts showing observed wi eking heights and calcul ated wettability, i.e. , cos 8 , 
are presented in Table 6-10. Because the test rig was designed to provide only an approxi- 
mate measurement of wicking height the values shown in Table 6-10, which are accurate 
to ±0. 64 cm, should be treated as a qualitative comparison between the test fluids. From 
the results it is clear that ethane has the same wetting characteristics as Freon-113 and 
acetone with respect to stainless steel screen. A wicking height of 1. 5 in. (3. 81 cm) was 
observed for ethane, which indicates a wettability factor of about 0. 5, comparable to that 
of Frepn-113 and acetone. 

It was observed during testing with all three fluids that the insertion of part of a dry 
sample of a single layer of screen into the fluid did not result in the fluid rising up along the 
wick. Rather it was found that the sample first had to be wet along its length then held’ 
vertically with one end in the fluid. In this way, drying would begin at the top of the wick 
and proceed to a certain point down the wick below which the wick continually remained wet. 

The situation in the diode heat pipe is quite different since there are many layers, 
of screen that make up the artery. The gaps formed between these layers promote capillary 
action that eventually wets the entire wick. Open air bench tests with a multiple layer screen 
wick dipped in acetone confirmed that the wick will draw fluid up along its length. This self- 
wetting property was also demonstrated with the multiple layer screen wick and ethane fluid 
using the test apparatus. 
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Table 6-8 Wick Sample Characteristics 


TYPE 

100 MESH 
STAINLESS STEEL 

CONSTRUCTION 

SINGLE LAYER 

1 BY 4 INCH 

(2.54 cm BY 10,16 cm) 

WIDTH OF 

0.0055 (0.1397) 

OPENING, INCH (mm) 


r, PORE RADIUS, 

0.00275 (0.0699) 

INCH (mm) 

' . • 

CLEANING 

PASSIVATED 

0787-047 ' 



Table 6-9 Fluid Properties 


. FLUID 

FREON-1 

ACETONE 

■ • ETHANE 

TEMP °F (K) 

86° (303 K) 

70°F (294 K) 

-120°F (189 K) 

. . DENSITY, LB/FT 3 

97 

48.9 

34 

: SURFACE TENSION 
LBF/FT 

0787-048 

' 0.0013 

0.00155 

‘ 0.00106 


Table 6-10 Wettability Test Results 



TEMP 
°F (K) 

WICK SAMPLE 

100 MESH 
PASSIVATED 

70 (294) 

0.75 (1.91) 

70 (294) 

1.75 (4.45) 

-123 (187) 

1.5 (381) 

70 (294) 

0.53 

70 (294) 

0.53 

-123 (187) 

0.46 






























Section 7 


CONCLUSIONS 

The development of spiral artery cryogenic thermal diodes has been advanced to in- 
clude the following accomplishments: 

• A liquid blockage (blocking orifice) design was built and qualified for the HEPP 
spaceflight experiment 

e A liquid trap engineering model was built and tested. Throughput (w-m) and for- 
ward mode conductance was similar to the liquid blockage design but reverse mode 
conductance and shutdown energy were higher 

• An analytical and experimental investigation of the startup dynamics of steel/ethane 
heat pipes with spiral artery wicks revealed that wick permeability is an important 
influence on the maximum heat load permitted during startup. Higher initial heat 
load capability (priming load) can be obtained by increasing the wick permeability. 

A modified composite artery wick should be built to evaluate design improvements 
aimed at increasing the priming heat load capability. It would include a coarse interior 
screen with relatively large and uniform gaps (for increased permeability), surrounded by 
a closure wrap of very fine pore screening (for maximum capillary pressure rise). 
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Appendix A 


HEPP Liquid Blockage Thermal Diode, Computer Model 


A-i 




0707-050 


200 


-Hh 


NODE 


DESCRIPTION 


*1 



202 


1 EVAPORATOR BLOCK 

2 EVAPORATOR WALL 

3 HP VAPOR (BOUNDARY NODE) 

4 CONDENSER WALL 

5 RESERVOIR 

6 RESERVOIR RADIATOR, ROOT 

7 MAIN RADIATOR (L 3 + L 4 ), ROOT 

8 MAIN RADIATOR (L 1 + L 2 ), ROOT 

9 PCM 

200 SPACE 

201 INTERNAL STRUCTURE 

202 EXTERNAL SHIELD 





MASS CAPACITANCE 



THERMAL MASS j 

NODE NUMBER 

(BTU/°F) 

(WATT-HR/°C) 




1 


.Uo/U 

.0459 

2 


.00557 

.00294 

4 


.0324 

.0171 

5 


.0160 

.00844 

6 


.0217 

■ .01 14 

7 


.078 

.0411 

8 


.174 

.918 

- 


.72 (SOLID) 

.38 

9 ■ 


108.6 (MELT) 

57.26 


- 

.956 (LIQUID) 

.504 
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CONDUCTION 




COUPLING VALUE j 

COUPLING NUMBER 

JOINING NODES 

(BTU/HR °F) 

(WATTS/°C) 

1 

1-201 

1 254 x 10 -4 

6.614 x 10" 5 

2 

3-201 

4.854 x 10" 4 

2.560 x 1 0' 4 

3 

1-2 

20.480 

10.801 

4 

2-3 

27.4 

14.45 

s 

3-4 

311,0 

164.02 

6 

4-5 

.773 

.408 

7 

4-8 

46.9 

24.73 

8 

7-8 

1.44 

.759 

9 

9-7 

10.000 

5.274 

10 

5-6 

2.75 

1.45 

RADIATION 



(BTU/HR °R 4 | 

(WATTS/°K 4 ) 

1 

1-201 \ 

2.242 xIO" 1 2 

689 x 1 0' 1 2 

2 

3-201 

6.18 x 10' 12 

1.9 x 10' 12 

3 . 

4-201 

3.259 x 10" 12 

10.02 x 10~ 12 

4 

5-201 

1.051 xIO' 12 

3.233 x IQ 42 

5 • 

8-201 

18.6 x 10’ 12 

578 x 10" 12 

6 

9-201 

5.044 x 10‘ 12 

15.51 x IQ' 12 

7 

6-201 

1.607X 10' 12 

4843x1 O' 1 2 

8 

8-202 

9.62 x 10~ 12 

2.96 x 10" 1 2 

9 

8-200 

9.52 xIO' 10 

29.3 x 10‘ 1D 

10 

7-200 

-10 

6.86 x 10 

21.1 x 10' 10 

11 

7-202 

6.93 x 1 0' 1 2 

21.3 x 10' 12 

12 

6-200 

94.1 x 10 42 

289 x 10' 10 

13 

6-202 

-12 

.95 x 10 ^ 

282 x 10' 12 
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INPUT VARIABLE 



T(255) 


T(256) 



RESERVOIR RADIATOR FIN LENGTH (IN), ... 

RESERVOIR RADIATOR SURFACE AREA (IN 2 ) 

RESERVOIR RADIATOR FIN THICKNESS (IN) 

PCM RADIATOR AR EA'(FOR L 4 >(IN?) ; • 

PCM RADIATOR FIN LENGTH (Lg) (IN) 

PCM RADIATOR FIN THICKNESS (IN) 

PCM RADIATOR AREA (FOR Lg) (IN 2 ) 

DIODE RADIATOR AREA (FOR Lg) (IN 2 ) 

DIODE RADIATOR FIN LENGTH (Cg) (IN) 

DIODE RADIATOR THICKNESS (IN) 

DIODE RADIATOR FIN LENGTH (FOR L, ) (IN) 

DIODE RADIATOR AREA (FOR Lj) (IN 2 ) 

PCM RADIATOR FIN LENGTH (L 4 ).(IN) 

— WIDTHOFPCM AND DIODE RADIATORS (IN) 

RESERVOIR’ RADIATOR EMITTANCE 
PCM AND DIODE RAD1ATQR EMITJANCE . 

T1, LOWER LIMIT OF MELT RANGE (°F) 

T2, UPPER LIMIT OF MELT RANGE (°F) 

LATENT HEAT OF PCM (BTU/LBM) 

MASS OF PCM MATERIAL (LBM) 

CP (SOLID) (BTU/LB °F) 

CP (LATENT) 

RATI = cP(SOLID) 

CP (LIQUID) 

RAT2 = CP (LATENT) 

CHECK VALUE TO DETERMINE WHAT CP TO USE (CPMCK) (BTU/°F) 

FRACTION OF PCM REMAINING 

WATT-HOURS REMAINING 

CAPO) - INITIAL VALUE 

COUNTER (MUST INITIALIZE TO 0.0) 

OPTION CODE FOR HEATER ROUTINES 
f +1 ACCEPTS HEATER CYCLING ROUTINE 
0 REJECTS BOTH OPTIONS 
[ -1 ACCEPTS DUTY CYCLE ROUTINE 
% DUTY CYCLE ON PCM 
TIME SPAN FOR PCM (SEC) 

% DUTY CYCLE ON RESERVOIR 
TIME SPAN FOR RESERVOIR (SEC) 
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OUTPUT VARIABLE 


.-DEFINITION 


T1205) - ' ' 

HEAT INPUTTO THE PCM.(BTU/HR) 

T(206} 

HEAT INPUT TO THE RESERVOIR (BTU/HR) 

T1208J 

NET HEAT AT DIODE PIPE (BTU/HR) 

T(209) 

CURRENT MASS IN RESERVOIR (GMS) 

T(210) 

MASS IN BLOCKED SECTION (GMS) 

T(211) 

EVAPORATOR LENGTH (INCHES) - 

T(212) 

BLOCKED LENGTH (INCHES) 

T(214) 

CONDUCTANCE BETWEEN EVAPORATOR & CONDENSER (BTU/HR°F) 

' T(215) 

LATENT HEAT OF WORKING FLUID (BTU/LBM) 


LIQUID DENSITY OF WORKING FLUID (LBM/FT 3 ) 

T{217) 

S MASS IN RESERVOIR (GMS) 

T(220) 

EVAPORATOR TEMPERATURE RATE OF CHANGE (°F/HR) 

T(221i 

P.CM TEMPERATURE RATE OF CHANGE (°F/HR) 

T(222) 

CONDENSER TEMPERATURE RATE OF CHANGE (°F/HR) • 

T(223! 

RESERVOIR TEMPERATURE RATE OF CHANGE (°F/HR) 

T(226) 

REVERSE MODE ENERGY (WATT-HRS) 

T(227) 

% RESERVOIR CHARGE REMAINING 

T(228) 

OUTPUT INTERVAL (HRS) 

T(229) 

PREVIOUS PRINTOUT TIME (SEC) 

T(230) 

FIN EFFECTIVENESS FOR RESERVOIR RADIATOR 

T(235) 

FIN EFFECTIVENESS FOR L 3 PORTION OF PCM RADIATOR' 

T(240) 

FIN EFFECTIVENESS FOR DIODE RADIATOR FOR 1^ 

T(243) 

FIN EFFECTIVENESS FOR DIODE RADIATOR FOR L, 
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o n 


F ILEO T15CLC1 FORTRAN PI 


GRUMMAN 


DATA 


SYSTEM 


SUBROUTINE CLC1 

THIS SUBRQ.VTANE..I$ -CABLED AT.. THE START . QF_ EVER T_ JLTER A TIQ.N . . 

COMMON ST G( 5000 | • NP0S(50), L 1 < 50 ) . L2(50>. NLAST(50), FCT(50) 

' COMMON KILL « NTABS. NMTR IC .METPR . I HEAT .SIG.XJ1 «X J2.X J3 
COMMON BUFF (20) .T( 1000) .CAP ( 1000) »COND(2000) .RAO (2000 ) .Q( 100) . 

1 NCONO( 2.2000 ) , NRAD ( 2.2000) , NI N{ 1 00 0 I . NBN ( 250 ) , KE Y ( 5 ) . NQ ( 100 > . 

2 KAY( A) ,FLUX( 1000) 

COMMON TIME. OTIME, FT I ME . OUM, 

1 KNIN. KNBN. N4» N5. KOOE, KNQ. 

2 KOUT. KOUNT, NOSRT 

REVISED TIROS-N CRYO EXPERIMENT NETWORK 

IF(TIME.GT.O.O) GO TO 132. 

SUMQ=1.0 . ' 

BKL=0.0 
HPCOND=25.ia 
XLAT=1 .0 

RHOL= 1.0 • __ 

AK=Q . 0 


TO 133 


TO 134 


TTA0001 0 
TTA.00020 
TTA00030 
TTA00040 
TTA00050 
TTA00060 
TTA00070 
TTA00080 
TTA00090 
TTAOOIOO 
TTAO 0110 
TTA00120 
TTAO 0130 
TT AO 0140 
TTA00150 
TTAO0160 
-T-T AO 01-70- 
TT AO 0180 
TTA0019Q 
TTAO 0200 
TTA00210 
TTA00220 
TTAO 0230 
TTA00240 
TTA00250 
TTAO 0260 
TTA00270 


132 AX=TIME-AK* 1.667 
IF( AX .LT. I . 667). GO 
AK— AK4- 1.0 

133 IF(AX.LE. 0.56335) GO 
Q(4)=0.0 
Q(5>=0.0 
GO TO 135 

134 Q(4)=,417*T(232)*SIN(5.5738*AX) 

Q{5)=.417*( T( 239)4T(24S) ) *S I N< 5 .57384AX ) 

135 T ( 230 ) =T ANH ( 2.390E— 06*T( 231 )*(T(248)*< (T( 6 ) +460 . »**3 ) /T( 233 ) ) **.5) TTA00280 

A/(2.390E-06*T(231 »*(T(248»*( ( T ( 6 ) +460 . 1**3 > /T( 233 ) >**.5> TTA00290 

RADI 12>=1.061E-1 1*( T< 232>*T( 2^0) ) * TTA00300 

RAD{ 13 ) - 1 .0 71E— l 3* ( T ( 232 >*T ( 230 ) ) TTA00310 

T ( 240 ) =T ANH ( 2.390E~06*T( 241)*<T( 249) *( CT< 8) +460 . >**3 ) /T{ 242) > ** .5 ) TTA00320 
C/( 2. 390E-06*T{ 241 )*(T( 24 9)*((T( 8)4-460. )**3 ) /TC 2421 )$*.S) TTA00330 

T(243)=TANH(2.390E-06*T(244)*(T(249)*( (T( 8 ) +460 . ) *43 ) /T (242)1*4.5) TTAO 0340 
D/(2.390E-06*T(244)*(T(249)*((T(8>*460. ) **3 ) /T ( 242 ) ) ** . 5 ) TTA00350 

RAD(8)=9.522E~14*{ T(239)*T{ 240 > +T ( 243 ) *T< 245 > > TTAO 0360 

RAD(9)=9.427E-12*( T(239)*T(240)4-T(243> *T(245>) TTA00370 

RAD(5) = 1.785E-13*( T ( 245 ) 4-T( 239) ) TTA00380 

RAD( 7) =1 ,73SE-13*T( 232) TTA00390 

CAP(6)=.0214*(T(232)*T(233) ) TTA00400 

CAP (8)=.0H1 44 ( ( T { 245 ) 4-T ( 239) ) *T{ 242 >4-2 .906 > TTA00410 

T( 3> = ( COND(4»*T( 2) 4-COND( 5) 4T(4 ) 4-C0ND( 2 ) *T( 20 1 >+RAD(2 ) *i ( T ( 20 1 ) 4-460TTA0 C420 
A. ) <=*4) ) /( COND ( 4 ) >COND ( 5 ) 4-COND<2 ) + (RAD(2)A( ( T(201 ) 4-460 . )**3) ) ) 

IF(T{264))401 .400.402 


401 


403 

404 

405 


CYCLE=TIME*3600./T( 272) 
ICYCLS— CYCLE 
CYTI ME=CYCLE- IC YCLE 
IF(CYTIME— T (271 ) >403.403.404 
Q( 1 )=T ( 269) 

GO TO 405 

Q( 1)=0.0 

RYCLE=TIME43600./T(274) 
IRYCLE-RYCLE 
RYTIME=R YCLE- IR YCLE 
IF(RYTIME-T<273) >406.406,407 


TTA00430 
TTAO 0440 
TTA00450 
TTA00460 
TTA00470 
TT AO 0480 
TTA00490 
TTA00500 
TTA00510 
TTAO 0520 
TTAO 0530 
TTA00540 
TTAO 0550 
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FILEO T15CLC1 FORTRAN PI 


GRUMMAN 


DATA 


SYSTEMS 


406 a(Z)=T(270> 

GO 10.-400 

407 Q(2)=0.0 
GO TO 400 

402 IF(T(4). Lt » T( 26 S) ) Q(1)=T(269) 

IF(T(4).GE.T(265).AND.T(4) .LE.(T( 265)+T<266 ) ) ). Q(1>=Q(1> 
. IF(T(4>.GT.(T(265)+T(266)>> Qtl>=0.0 ' ■ 

IF(T{5) .LT_»JTJL267> ) Q<2» = TC27.0) 

IF ( T( 5 1 «GE • T( 267 ) . AND .T < 5) .LE. ( T 1 267> +.TC 268 ) ) ) Q(2)=Q(2> 
IF ( T( 5 1 »GT » (T (267 ) 4TI268 ) I ) Q(2»=0o0 
400 T(205)=Q<1) 

T(206)=Q(2) 

RESI=1 .376/453.5924 
EVPL=4./12. 

IF(BKL.GT.O.O) GO TO 600 

C0ND(4)=27.4 

COND< 51=311.0 

HPCOND=l./( ( 1 ./COND (4) 1 + <1. ✓CONDI 5))) 

600 IF<TIME.GT.O.» GO TO ISO 
T20=T(2) 

T90=T(9> 

T40=T(4) 

T50=T{5> 

QREV=0. 

150 IF(T( 3). GT . T ( 2 ) ) GO TO 100 


BLKM=0. 

BKL=0. 

DR£M=0. 

RESM=R£S1 
GO TO 14 

• 100 IFCRE5M.GT.0. ) GO TO 200 
DREM=0 . 

RESM=0.0 

BLKM=RESl 




IF(T(4) .LT.TC2) ) GO TO 14 


GO TO 8 


200 T5R=460.+T(S) 

T20IR=460.«-T<201) 

SUMa=Q{4)-RAD(4>*(T5R**4.-T201R4*4.>-COND{ l 0 )*{ T( 5 ) -T< 6) > 
TT=T(5» 

XLAT=PR0P(6.6»TT) 

RHOL=PKOP( 1 .6.TT1 

DREM={ SUMQ4DT IME-CAP ( 5 I’M T ( 5 )-T50 ) ) /XL AT 
RESM=RESM-DREM 
BLKM=8LKM40REM 
8 IF ( BLKMoLT .0.0) BLKM=0.0 
BKL=BLKM/( RHQI.4 1 . 2566E-4 ) 

IF(BKL .LT.O .0 ) BKL=0.0 
IF(BKL.LT.EVPL) GO TO 300 
BKLT = BKL-EVP1_ 

I F( T ( 4 ) .LE .T( 2) ) . GO TO 14 
COND (4 1=9.1556-04 
CON0(5)=174.0 

QREV=9 » 155E— 044 (T(41— T(2))+QREV 
GO TO 14 


TT AO 0560 
TTAQ057 0 
TTA00580 
TTA00590 
T-TA0 0600 
TTA0061 0 
TTA00620 
TTA00630 
TTA00640 
TTA00650 
TTA00660 
TTA00670 
TT AO 0660 
TTA00690 
TT AO 0 70 0 
TTA00710 
TTA00720 
TTA00730 
TTA00740 
TTA00750 
TTA00760 
TTA00770 
TTA00780 
TTA00790 
TTA00800 
TTA00810 
TTA00820 
TTA00830 
TTA0C840 
TTAOOC50 
TTA00860 
TTAOC870 
TTA00880 
■TTA00890 
TTA00900 
TTA0091 0 
TTA00920 
TTA00930 
TTA00940 
TTA0QS50 
TTA00960 
TTA00970 
TTA0C980 
TTA00990 
TT AO 100 0 
TTA01010 
TTA01020 
TTA01030 
TTA01040 
TTAOIOSO 
TTA01060 
TTA01C70 
TTAOIOSO 
TTA01090 
TTA01100 
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FILEO T15CLCI FORTRAN PI 


GRUMMAN DATA 


S Y S‘ T E M S 


300 

IF(T(4).LE.T(2)) GO TO IA 




. TTAO 1110 


*BKk + SflBJX »£ JL226 * BJSJ^JLAN.HI. LL2, ST*B<t ) 

“-T.TA.CI Jifl, 


COND<4)~HPCCN0 




TTAO 1130 


COND<5)=174.0 




TTAO 1140 


QREV=HPCOND*( T ( 4 )— T ( 2 ) M-QREV 




, TTAO 1150 

14 

CONTINUE 




TTAO 1 160 

PCM 

CALCULATIONS, T l=LOWER MELT TEMP. T2-=UPPER 

MELT 

TEMP. 

XLAT=LAT 

HT TTAO 1170 

SLOST^ WT PCM* CPSLD— CP SQUID. RATI— CP(MELT) 

/CP' 

SOLID 

. RAT2=CP(LQ) TTA01180 

/CP ( MELT 1 , CPMCK= CP CHECK VALUE, MCP SOdO. 

KCP(LQ) .LT. 

MCP(MEUT) 

• TTAO 1190 


COMMON /PCM/ R AT i,Tl,T2, CPMCK 




TTAO 1200 


CPMCK=T(2S7) 




TTAO 1210 


T1=T(250> 




TTAO 1220 


T2=Tt25l| 




TTAO 1230 


'RAT 1— T C 255 > 




TTAO 1240 


RAT2=T< 256) 




T.TA012S0 

. 

TMID=TI/2. ♦ T2/2. 




TTAO I 260 

T(9)=PCM NODE 




TTAO 1 270 


IF(TC9 > .LT. Tl ) GO TO 10 _ 




TTAO 1280 


I'F(T<9) .LT.TMID > GO TO 20 




TT AO 1290 


IF(T(9> .LT.TZ) GO TO 30 




TTA01300 

- 

IF<CAP<9) ,GT. CPMCK) CAP(9)=CAP|9)*RAT2 




TTA01310 


GO TO 98 




TTAO 1320 

10 

1F(CAP(9) ,GT. CPMCK } CAP <9 ) =CAP( 9) /RATI 




TT AO 1330 


GO TO 98 




TTAO 1340 

20 

IF ( CAP ( 9 J .LT. CPMCK) CAP ( 9 ) — CAP( 91 4RAT1 




TTAO l 350 

. 

GO TO 98 




TTAO 1360 

30 

IF(CAP{ 9 ) .LT. CPMCK ) C AP 1 9 > =C AP C 9 ) /RA T 2 




TTA01?70 

98 

CONTINUE 




TTAO 1380 

A 

T{20Q)=SUMQ 




TTAC1390 

• 

T { 209) =RESM4453. 5924 




TT AC 1400 


T( 21 0 ) =BLK M *4 53 • 59 2 A 




TTAO 1410 


T<211)=C0ND<4) 




" TTAO 14 20 


Tf 212)=BKL*12. 




TTAO 1430 


Tl21Al=CONOCS) 




TTAO 1440 


T(215)=XLAT 




TTAO 1450 


T( 216) =RHOL 




TTA01460 


T ( 21 7 ) =DREM*453 . 5924 




TTA01470 


T1220)=CT(2)-T20)/DTIM£ 




TTAO 1480 


T( 221 ) = ( T{ 9 1— T90 )/DTI ME 




TT AO 1490 


T ( 222 ) = <T(4)— T40) /OTI ME 




• TTAO 150 0 


T ( 223 ) = C T( 5 ) — T50 ) /DT I ME 




TTAO 151 0 


T ( 226 )=QREV*DTIME/3 .4 13 




' TTAO 1520 


T ( 227 ) =T (20 9) /(RES 14453. 5924) 




TT AO 1 530 


T ( 228) =KOUT 




TTAC 1540 


T ( 229) -TIME43600. 




TTAO 1550 


T20=T( 2) 




TTAO 1560 


T90=T (9) 




TTA0157C 


T40=T ( 4 > 




TTAO 1580 


T50=T(5) 




TTAO 1 590 


T(263)=CAP<9> 




TTAO 1 600 


T ( 234)— Q<4) /3.4X3 




TTA01610 


T(235>=a(5)/3.413 




TT AO 1620 


RETURN 




TTA01630 


END 




TT AO 1 640 
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FILEO TCASEl DATA PI GRUMMAN DATA SYSTEMS 

1. REVISED TRANSIENT MODEL FOR TIROS-N CRYO EXPERIMENT 

3. 

4 » BY BRAZE AT RESERVOIR SADDLE 
S. NOOE 3 ITERATED 


0 

0 


1 

0 



0 

0.0 


5 

.OOOE- 

os 

5. OOOE 

0 

I- 

9 

• OJlftE 01 

0 

0 

0 

2-9 

• 000E 

01 

0 

0 

0 

- 3- 

9 

*400E 

01 

0 

0 

0 

4- 

9 

• SOOE 

0 1 

0 

0 

0 

5- 

9 

• 800 E 

01 

0 

0 

0 

6-9 

.eooE 

0 1 

0 

0 

0 

8- 

9 

>B00E 

01 

0 

0 

0 

9- 

8 

.oooe 

0 1 

0 

0 

00 

0 

0 

.0 


0 

0 

0 

200- 

4 

. 600E 

02 

0 

0 

0 

201 

7 

.OOOE 

0 1 

0 

0 

0 

202- 

2 

• OOOE 

02 

0 

0 

0 

205 

0 

.0 


0 

0 

0 

206 

0 

.0 


0 

0 

0 

208 

0 

.0 


0 

0 

0 

209 

0 

.0 


0 

0 

0 

210 

0 

.0 


0 

0 

0 

21 1 

0 

.0 


0 

0 

0 

212 

0 

.0 


0 

0 

0 

214 

0 

.0 


0 

0 

0 

21S 

0 

.0 


0 

0 

0 

216 

0 

.0 


0 

0 

0 

217 

0 

.0 


0 

0 

0 

220 

0 

.0 


0 

0 

0 

221 

0 

.0 


0 

0 

0 

222 

0 

.0 


0 

0 

0 

223 

0 

.0 


0 

0 

0 

226 

0 

.0 


0 

0 

0 

227 

0 

.0 


0 

0 

0 

228 

0 

.0 


0 

0 

0 

229 

0 

.0 


0 

0 

0 

230 

8 

.OOOE— 0 1 

0 

0 

0 

231 

4 

*690E 

00 

0 

0 

0 

232 

22.00E 

00 

0 

0 

0 

233 

2 

.OOOE— 02 

0 

0 

0 

234 

0 

.OOOE 

00 

0 

0 

0 

235 

0 

.oooe 

00 

0 

0 

0 

239 

15.00E 

01 

0 

0 

0 

240 

8 

.OOOE— 0 1 

0 

0 

0 

241 

8 

. 660E 

00 

0 

0 

0 

242 

2 

• OOOE- 

-02 

0 

0 

0 

243 

8 

• OOOE-O 1 

0 

0 

0 

244 

8 

.660E 

00 

0 

0 

0 

245 

15.00E 

01 

0 

0 

0 

248 

9 

.OOOE-O 1 

0 

0 

0 

249 

9 

.OOOE-O 1 

0 

0 

0 

250- 

-1 

♦ 315E 

02 

0 

0 

0 

251- 

1 

.305E 

02 

0 

0 
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FILEO TCASE 1 


DATA. . 


DATA 


SYSTEMS 


PI 


G R U -M M A N 


0 

252 

6 

»030E 

01 

0 

0 


0 

253 

3 

.OOOE 

.09 

0 

0 


0 

254 

4 

• OOOE-O 1 

0 

0 


0 

255 

1 

. 508E 

02 

0 

0 


0 

256 

8 

.8006-03 

0 

0 


0 

257 

1 

• OOOE 

02 

0 

0 


0 

260. 

1 

.OOOE 

oo 

0 

0 


0 

261 

1 

.OOOE 

00 

.0 

0 


0 

262 

0 

.0 


0 

0 


0 

263 

0 

.0 


0 

0 


0 

264-1 

.0 


0 

o 


0 

265-9 

• 400E 

0 1 

0 

0 


0 

266 

1 

» 800E 

01 

0 

o • 


0 

267-8 

• 500E 

01 

0 

0 


0 

268 

1 

* 800E 

01 

0 

0 


-0 269- 

*0 

TOOOE - 0 2 

0 

0 


0 

270 

0 

.OOOE 

00 

0 

0 


0 

271 

4 

.OOOE 

-01 

0 

0 


0 

27 2 

1 

.OOOE 

01 

0 

0 


0 

273 

1 

.OOOE 

00 

0 

0 


0 

274 

1 

.OOOE 

01. 

0 

0 


1 1 100 

0 

0 

.0 


0 

0 


0 

1 

8 

.7 OOE 

-0 2 

0 

0 


0 

2 

5 

. 570 E 

-03 

0 

0 


0 

3 

6 

. 150E 

-02 

0 

0 


0 

4 

3 

. 240E 

-02 

0 

0 


0 

S 

5 

.OOOE 

-0 3 

0 

0 


0 

6 

1 

• OOOE 

00 

0 

0 


0 

8 

1 

.OOOE 

00 

0 

0 


0 

9 

1 

. 593E 

00 

0 

0 


11100 

0 

0 

.0 


0 

0 


0 

1 


1 

201 

1.254E- 

04 

0 

0 

2 


3 

20 1 

4.854E- 

04 

0 

0 

3 


1 

2 

2.048E 

01 

0 

0 

4 


2 

3 

27.40E 

00 

0 

0 

5 


3 

4 

31.10E 

01 

0 

0 

6 


4 

5 

7.730E- 

01 

0 

0 

7 


4 

8. 4.690E 

01 

0 

0 

10 


5 

6 

2.750E 

00 

0 

0 

11 


8 

201 

0.0 


0 

0 

13 


6 

201 

0.0 


0 

11100 

0 


^0 

0 

0.0 


0 

0 

1 


1 

20 1 

2.242E- 

12 

0 

0 

.2 


3 

201 

6.180E- 

12 

0 

0 

3 


4 

20 1 

3.259E- 

12 

0 

0 

4 


5 

201 

1.05 1E- 

12 

0 

0 

5 


8 

201 

1 .OOOE— 

09 

0 

0 

6 


9 

201 

S.044E- 

12 

0 

0 

7 


6 

201 

l.OOOE- 

09 

0 

0 

8 


8 

202 

l.OOOE- 

09 

0 

0 

9 


a 

20 0 

l.OOOE- 

09 

0 

0 

12 


6 

200 

1 .OOOE— 

09 

0 

0 

13 


6 

202 

1 .OOOE— 

09 

0 

11100 

0 

< 

0 

0 

0.0 


0 

0 

1 


9 0 .0 OOE 02 




0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 
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RHPKOD (JC1BILITY OF THE 
ORIGINAL PAGE IS POOR 


FILEO TCASEl DATA 


PI 


0 

0 

0 

0 

11100 

X HOC 
0 

1 . 

2 . 

3 . 

A. 

5. 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

moo 

ll-ioo 

11100 

moo 

moo 

moo 

moo 


2 

3 

A 

5 

0 

0 


5 O.OOOE 00 
• --1- 0 .QQCE Of) 

6 0.0 
8 0.0 
0 0.0 
0 


0.0 

1 

2 

4 

5 

6 
8 
9 
3 

200 - 

201 

202 - 

262 

0 

0 

0 

0 


2000 0 

5. 0006-03 2.40OE 01 


7.7006 01 
7.700E 01 
7.7 OOE 0 1 
7.7006 01 
7.7006 01 
7.700E 01 

7.7006 01 

7.7006 01 
•4,6006 02 

7.0006 01 

2.0006 02 

0,0 

0,0 

0,0 


C 

a 


o.o 


Oi 

o. 


o 

o 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 


0 

0 


GRUMMAN DATA 


SYSTEMS 
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Appendix B 


Subroutines CLC1 and CLC2 for T15 Computer Code (Priming Analysis) 


B-i 



ssssvs- 


RELEASE 2.0 CLC1 DATE = 79075 16/46/24 


SUBROUTINE CLC 1 

C THIS SUBROUTINE IS CALLED AT THE STAHT OF EVERY ITERATION 

COMMON STG(SQOO), NP0S(5 0), Li (50) , 12(50), ill AST (50), FCI(50) 

COMMON KILI, NTABS,NHTRIC,l!ETP8„IBEAT,SlG,X01,XJ2»XJ3 

COMMON BUFF (20) , I (1000) , C&£> (1000) ,COND(2000),BAD (2000) ,Q(100) , 

1 NCOND{2, 2000), NBAD (2,2000) ,NIN(1 000) , NBN (250) , KEY (5) ,NQ(10G) , 

2 ' KAY(4) ,FIUX ( 1 00 C) , DISl (2000) ,DIS2{2000) ,ISC1 (2000) ,ISC2 (2000) 

COMMON TIME, DTIME, FIX ME, DUH, 

1 KNXN, KNBH, N4 , N5, KODE, KNQ, 

2 KOUT, KO.UNT, NOSfiT 

DIMENSION A(43) ,TB(8) ,P{10) ,RHOV{10) ,XMUV (10) ,HFG(10) 

C CCND (54) =PV, SAT. PEESSUBE FOR T<24) , PSF 

C COND (55)=DE1PP, PRESSURE DIFF- AVAILABLE FOR GAP PRIMING, PSF 

C ' CCND (56) =BEOl, LIQUID DENSITY OF NODE 28, PCF 
C CCND (57)=PGAP, GAP ERESSORE,PSF 

C COND (58) =PCGAP, CAPILLARY PRESSURE OF GAP, PSF 

C COND (59) = KVAP, FLUID IDENTIFICATION FOR PROP SUBROUTINE 

C 'COND (60)- B, EXPONENT FOB MLI PARALLEL CONDUCTANCE CALC 

C COND (61)= CPV, VAPOR SPECIFIC HEAT 

C CCND (62)= CPI, LIQUID SPECIFIC HEAT 

C COND (63)= TGAP, GAP THICKNESS, INCHES 

C COND (64)= HGAP, GAP HEIGHT ABOVE LOH POINT OF CONDENSER 

C COND (65) -COND ( 107).= INITIAL VALUES OF ConD 1-43 RESPECTIVELY 

C COND ( 1 08) =1-01 SKIPS INITIAL BOOKKEEPING 

C CCND (109)=T(24) 

C COND (110) =T (32) 

C COND (118) =DPI, LIQUID PRESSURE DROP, PSF 

C CCND (1 19)=DABT, ARTERY O.D-, INCHES 

C COND (120) =APOR, POROUS AREA, SQ- IN. 

C CCND (121)=DPOK, POROUS MAIL. HYD. DIA-, INCHES 

C COND (122) =DI, INSIDE DIA. OF PIPE, INCHES 

C K VAP= 6 FOB ETHANE 

C COND (123) = XMLI, COB- FACTOR FOR MLI RADIAL RAD. LINKS 

C COND (124) =C0ERECTI0N FACTOR (FRACTION) TO ARTERY PER8EABILITY*AHEA 

C COND (125) =ARTSR¥ PE 8ME ABILITY* AREA, FX**4 

C COND (126) =CO ERECTION FACTOR FOR Mil PARALLEL CONDUCTANCE 

C COND (127) =TOTAL HEAT DECREMENT IN HA TIS To FIND PRIMING LIMIT 

N1=COND(108) 

IF (N 1- FQ. 1) .GO TO 2 
K?AP= COND (59) 

B = CCND (60) 

CPV= CCND (61) 

CPL= COND (62) 

TGAP=C0ND(63)/12. 

HGAP=COND (64)/ 12. . 

C STORE INITIAL CONDUCTANCE VALUES FOR ALL REGULAR CONDUCTORS 
DO 1 1=1,43 
1 COND (1+64) =CCND (I) 


T1500010 
T 1500020 
T 1500030 
T 1500040 
T 1500050 
T1500060 
T1500070 
T 1500080 
T 1500090 
T 1500 100 
T 15001 1 0 
T 1500120 
T150013Q 
T 1500140 
T 1500150 
T 1500 160 
T 1500170 
T 1500130 
T 1500190 
T 1500200 
T 15002 10 
T 1500220 
T 1500230 
T 1500240 
T 1500250 
T 1500260 
T1500270 
I 1500280 
T 1500290 
T 1500300 
1 1500310 
T 1500320 
X 1500330 
T 1500340 
T 1500350 
11500360 
T 1500370 
X1500380 
11500390 
T 1500400 
11500410 
T 1500420 
T 1500430 
T1500440 
T 1500450 
T 1500460 
T 1500470 
T 1500480 
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REIEAS.E 2.0 


C1C1 


DATE = 79075 


16/46/24 


COND { 108)=1- 01 

G "^ITER VAEOR PRESSURE DROP Bl FACTOR INPUT AS COND flit) 
PFAC=C0ND{11 1) 

DO 7 1=38,43 
7 GO NO {I) = CO ND (I) /PE AC 

C- ALTEB MLI EFF EHITTA NCE B X FACTOR INPUT AS COND (123) 
~XHII=CCND £723) 

DC 9 1=1,5 

HAD (I) =XMLI*RAD{I) . 

'• % RAD (I* 10) =XM1I*RAD (1+10) 

C- RECOMPUTE CONDUCTORS 1 THRU a 

2 EG 5 1=1,43 
5 A {I) =CCND (1+64) 

DO 12 1=1,8 


T 1500490 
T1500500 
T 1500510 
T 1500520' 
T 1500530 
T 1500540 
T1500550 
T 1500560 
T1500570 
T 1500580 
T 1500590 
T 1500600 
T150061 0 
— T-15 00620 - 


1 2— A-{-I) =CO ND (1 2 6) *A (X) 

DO 3 1=1,4 

TB(I) = <{ (T{I) *459.6) **4+{T (1+ 1) +459.6)**4) /2.) **.25 
3 COND{I)= A{I)*{TB{l)/36Q„) **B 
DC 4 3=5,8' 

TB {!) = ■(( {1(1+1) *459.6) **4+ (T{I+2) +459.6) **4)/2.) **-25 
_ - 4 COND {!) = A (I)*{TB [I;)/36Q_) .**B 
G gECOHEUTE FLUID CONDUCTORS 38 THRU 53 
TV=T{24) 

BFGY=BROP{ 6 , KVAP,TV) 

C TEMP TV MOST BE IN DECREES F FCB PROP SUBROUTINE 
C2=CPI/HFGV 
DC 6 1=1,7 
J=I+ 1 7 

P(I) = EROP{ 8 ,KVAP„I{J)) 

HFG (I)=PROP { 6 , KVAP, 3 (a) ) 

RHOV {I) = EBGP(3„KVAP,I( j)) 

6 XHUY {I) = PRO.R {4 ,K7AP,T{ J) ) 

DO 8 1=1,6 
J=I+37 

HFGB= { BFG (I) -fHFG (It 1) ) /2. 

RHOV E= {RHOV (I) +BHOV {1+ 1) ) /2. 

XMUVB= (XMU V (I) +XMO? (1+ 1) ) /2- . 

DPVDT=HHOVE*HFGB/ {T (1+ 17) +459. 6) *778. 

8 COND { J)= A (J) *BEVpT*KFGB*RHOVB/XMUVB 
COND (53) =C2*CQND{17)*{3 (11)-T (18)) 

COND {52)=C2*CQND(i8)*(3{12)-T (19)) 

COND (51) =C2* {COND {2 4) * (1(2 5) -T { 1 8 ) ) +COND {3 1) * {T(25|-T(32} ) ) +COWD 
1 (S3) 

CCND{50) =C2* {C-OND (25) * (T{26)-T{19)) +CQND{32) * (T { 26) ~T (32) ) )+COHD 
1 (52) +COND (51) 

COND {49) =C2* {COND {2 6 ) * {T {27) - T (20) ) +COHD {33) * (T (27) -T (32) ) ) +COND 
1{50) ‘ 

CCND (48) =C2* (GQND (2 7) * {T{28)~T {21) ) *CQfjD ( 34 ) * {T (28) -T (32) ) ) +COND 


T1500630 
T 1500640 
T 1500650 
T1500660 
T 1500670 
T 1500680 
T1500690 
T15007 00 
T 1500710 
T 1500720 
T1500730 
T 1500740 
T 1500750 
T 1500760 
T1500770 
T 1500780 
T 1500790 
T 1500800 
T15G0810 
T 1500820 
T 1500830 
T 1500840 
T 1500850 
T 1500860 
T1500870 
T1500880 
T 1500890 
T1500900 
T150G910 
T1500920 
T 1500930 
T1500940 
T 1500950 
T1500960 


B-2 




JEIEA3E 2.0 , CLC 1 DATE = 79075 16/46/24 


1 (49) 

COND (47) =C2* (CCND (2 8) * (T (29) — T (22) ) +CCND (35) * (T (29) -T (32) ) ) +COND 
1 (48) 

COND (45)=C2* (COND (22) * (T{23)-T (1 6} ) +COND (23) * (T (24) -T (17) ) ) /2- 
CCND (44) =CCND ( 45) 

COND (4fi)=C2* (COND (2 9) * (T(30) -I (23) ) *COND (36) * (T(30) -T(32) ) ) +CON D 
1 (47) -CCND (45) 

PY=PHCP(8,KVAP,TV) 

RH0L=PH0P(1,KVAP,T(31) ) 

STEN=EROP (5, Ky AP,T(32) ) 

DART=C0ND(119) 

E2=C0NE ( 122) 

NVAP= (DI-DAET) /24. 

STENV=EROP (5,KVAP,TV) 

PL=PV-BHOL*HGAP-2.*£XENV/HVAP 

ECgAp=ST2n*2./TGAP 

RH7- PBCP (3,KVAP,T (2 4) ) 

HG=PRCP(6,KVAP,T(24J) 

DPVDI=BHV*HG*778./ (T (24) <-459. 6) 

PGAP=EV*DPVDT* (T (32J-T (24) ) 

COND (54) =PV 
CCND (E6)=BH0I, 

COND (57) =PGAP 

CCND (58) -=PCGAP 

D ELPP= PI * P CG AP-PG A P 

COND (£5)=DELEP 

IF (CO ND (43) . EQ. 0„ ) GC TC 11 

IF(DEIPP. GT. C.) GO 10 10 

DELH= — DELPP/RHOL*! 2_ 

WRITE (6, 100) DEXH 
IF(DElH.LE.HGaP) WRITE (6,102) 

XG1=24.*STEN/ (PGAP-FL) 

IG0=2 4.*STEN/(EGAP-EV) 

WRITE (6,103) TGI 
WBIT E (6, 104) TGO 
GO TO 11 

10 WHITE (6,101) 

11 CONTINUE 

100 FORMAT (' ABTERX WILI NCT PRIME IN 1-G, EXTRA HEAD REQ- = *,F9. 5," 
1 INCHES*) 

101 FORMAT (' ARTeBX PRICES IN 1-G*) 

102 FORMAT (' AETERT PRISES IN 0-G') 

103 FORMAT {* MAX GAP WIITH FOR 1-G PRIMING - « e F6.5,« INCHES*) 

104 FORMAT!' MAX GAP WIDTH FOR 0-G PRIMING = *,F6»5 # « INCHES') 

RETURN 

C DEBUG INIT, SUBCHK 

END 


T1500970 
T 1500980 
T1500990 
T 150 1000 
T 1501010 
T 1501020 
T1501030 
T 150 1040 
T1501050 
T1501060 
T1501070 
T 1501080 
T 1501090 
T1501100 
T 15 01 1 10 
T 150 1 120 
T1501130 
T1501140 
T 150 1150 
T15Q1160 
T1501 170 
T 1501 180 
T1501190 
T1501200 
T 150 12 1 0 
T 1501220 
T 150 1230 
T1501240 
T 1501250 
T1501260 
T 150 1270 
T1501280 
T I50l 2 90 
T1501300 
T 1501310 
T 1501320 
T 150 1330 
T1501340 
T1501350 
T 1501360 
T 15 01370 
T 1501380 
T 1501390 
T 1501400 
T1501410 
T 1501420 
T 1501430 
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BEIEASE 2.0 


CLC2 


DATE = 79078 


10/46/29 


C 

c 


c 


c 


j3UBB CUTIHE CIC2 

THIS SUBBOUTIHE IS CALIED AT THE END OF EVEBY 
ITEHATICH BEFOBE TIHE IS IHCHEMEHTED 

CCHMCB STG(5000), NPCS{50) , 11 (50), 12(50), HIAST(50), FCT{50) 

COHHOH KI1L,HTABS,HETS2C,HETPH,IHEAT ,SIG,X J1 , X-J2,XJ3 

COtlfJOB B0FF(20) ,T (1000) , CAP (1000) ,COND{2000) ,BAD(2000) ,Q(100) , 

1 HCCND (2,2000) , HEAD (2,2000) ,HIN(1000) ,HBB (250) ,KEY (5) ,NQ(1 00) , 

2 KA X (4) ,F1UX (100 0) ,DIS1 (2000) ,DIS2(2000) ,lSci (2000) ,ISC2 (2000) 

CCHHCS TIME, DTIME , FTIHE, DUB, 

1 KKIK, 8HBH, H4, N5, KODE, KHQ, 

2 ROOT, KOOHT , HOSBT 

EIHE ESION A (43) 

THIS SOEBCUTINE PBR3ITS PABTIA1 AETEBY PBIHING 
DE1P P=COHD (55) 

IF (DEIPP.1E. 0. ) GO TC 10 

CPOB=CCND (12 1) . 

_P,I=3— 1-4.1-5-9 € ' 

BABT=CCND (119) 

APOB=CCHD(12C) 

AGAP=P I*CABT4*2_/4. -APCS-. 75* PI* (.01 6) **2 
TGAP=CCND (63) 

C 1=AGAP*TGAP**2./24£832. 

C2=AECB*DPCB**2./ 16 58880. 

C1F= (C 1+C2) *12- 
CIF=CIE*COHD (124) 

CCND (125)=C1F/12. 

H= 37 

2 IF(CCMD(H) .IT. .0001) GO TO 1 
COHD (H)=0. 

COHD (H-7) =COND (H-7) *.0821 
GC TC 9 
1 H=H-1 

IE (B.IE. 31) GC TO 10 
GO TC 2 

11 Q(1)=Q(1) -COHD (127) 43.413/2. 

C (2) =C (1) 

11=37 

IF(Q{1}.1B.G.) GO TO 10 
DO 12 1=31,37 

12 CORD (I)=COND (1*64) 

QEVAP- (Q(1)+Q (2) J/3.413 
HBITE (6,204) CEVAP 

9 CAIL STEADY 
KVAP=COHD (59) 

HGAP=CCND(64)/12. 

TGAP=COHD (63J/12. 

PV=PBCE(8,K7AP,T(24)) 

BHOL=EBOP ( 1,KVAP,T (31) ) 


T 1500010 
T 1500020 
T 1500030 
T1500040 
T 1500050 
T 1500060 
T1500070 
T1500080 
T 1500090 
T 1500 100 
T1500110 
T 1500120 
T 1500 130 
T15Q0140 
T1500150L 
T 1500160 
T 15001 70 
T 1500180 
T1500190 
T 1500200 
T 1500210 
T 1500220 
T 1500230 
T 15 00240 
T 1500250 
T 1500260 
T 1500270 
1 1500280 
T 1500290 
T1500300 
T 15003 10 
T1500320 
T 1500330 
T1500340 
T 1500350 
T1500360 
T 1500370 
T 1500380 
T 1500390 
T 1500400 
T 1500410 
T 1500420 
1 1500430 
T 15 00440 
T 1500450 
T 1500460 
T 15 00470 
T 1500480 
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DATE = 79070 


10/46/2 9 


STEN=EBOP ( 5, KVAP ,1 (32) )' 

PG-AP = E BOP £ 8, KV AP, I (32) ) . 

BBOV=IBOP (3,K¥AP,T{24) ) 

BFG =PBO P { 6, KVAP , T ( 24) ) 

BPYDT=HBCV*HFG*778. / (T (24) +459.6) 

PGAP=E¥+DPVDT* (I £32)-T(24) } 

ECGAP=STEN*2«»/IGAP 

EPL=Q. 

IHDL^EHOP (2,KVAP,1{31) ) 

C3=X EUI/BHC1/CCBE £62) /CLF/4. 1 8E08 
J=39— 8 
DC 4 1=1, J 

4 EPL=C3*COND (45+1) *CCND (1 1 1 *1) +DPI 
CCND (118) = EPI 
HVAP=(COND(122)-DAHT}/24. 

SlEnv^EEOP (5,KVAP,T (24)) 

EI=PV-BH0I*HGAP-2.*STEHV/¥VAP-EP1 
DELPP=EL+PCG AP— PGAP 

COND (54) = PV 
COND (55) =DE1PP 
CGHD (56) =BBCI 
CGHD (57) =PGAP 
CCHD (58)=PCGAP 
COND ( 109) =T (24) 

CCHD <1 10)=T (32) 

CALL BITE 
H=N-6 

IP (DEIEP.IT. 0.) GO TC 3 
HBITE (6,200) E,DE1EP 
GC IC 1 

3 BHITE {£, 201} H,DSLPP 

DEPZG=PSPSS0BE DlFP POB PBIHIHG IN ZEBO-G 

EPPZG=BELPP+BHOI*HGAP 

IP (DPEZG.LT. 0. ) GO TC 5 

HBITE (6, 202) E r DPPZG 

GC TC 1 

5 HBITE (6,203) 8,DPPZG 

IF (EElPP.LT. 0.) GO TC 11 

200 FCBHA7 (* ASTEBI P BIBBS IB 1-G HITH NODE’ #13, • PILLED, PHESSOBB 
IE^jFS.S,’ PSF*) 

201 PCBMAT £• ABTEBY HILI HOT PBIHE IN 1-G HlTH NODE* ,13, * PILLED, 
IS QBE E IFF =»,F9.5,» ESP*) 

202 POBHAT £« ABTEBY PBIMES IB 0-G 8ITH NODE * ,13, • PILLED, PfiESSOBE 
1 E2FFEEEHCE = *,F9.5,* PSP») 

203 FOBBAT ('ABTEBI SILL HOT PBIHE IH 0-G »ITH BODE *,I3,* FILLED, 
1S0BE riPP = ',?9.5,* PSP') 

204 POBHAI {//, * EYAPOBATOB POB IB CHANGED TO*,F9.5,» BATTS*) 

10 IF (DBIEP.LT. 0- .AND.C(I) .GT.O-) GO TO 11 


T 1500490 
T1500500 
T 15005 10 
T 1500520 
T 1500530 
T1500540 
T 1500550 
T 1500560 
T 1500570 
T1500580 
T 1500590 
T 1500600 
T1500610 
T 1500620 
T1500630 
T 1500640 
T1500650 
T 1500660 
T1500670 
T1500680 
T1500690 
T 1500700 
T150071Q 
T 1500720 
-T 1500730 
T 1500740 
T 1500750 
T 1500760 
T 1500770 
T1500*’80 
T 1500790 
T1500800 
T 1500810 
T 1500820 
I 1500830 
T 1500840 
T1 500850 
T 1500860 
DIPT 1500870 
T1500880 
PBEST1500890 
T1500900 
T 1500910 
11500920 
PBEST1500930 
T 1500940 
T 1500950 
T 1500960 
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OATS = 79078 


10/46/29 


CAXL BJTE 

8ETOBH 

EHD 


T 150097-0 
T1500980 
T 1500990 






